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Preface 


My  purpose  of  study  was  to  examine  the  non-rotating 
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Abstract 


The  assumption  of  a  non-rotating  planet,  common  in  most 
analytical  entry  trajectory  analyses,  has  been  shown  to 
produce  significant  errors  in  some  solutions  for  the  lifting 
atmospheric  entry  of  Earth.  This  thesis  presents  an 
investigation  of  the  validity  of  the  non-rotating  planet 
assumption  for  general  three-dimensional  Earth  atmospheric 
entry . 

In  this  effort,  the  three-dimensional  equations  of 
motion  for  lifting  atmospheric  entry  are  expanded  to  include 
a  rotating  planet  model.  A  strictly  exponential  atmosphere, 
rotating  at  the  same  rate  as  the  planet,  is  assumed  with 
density  as  a  function  of  radial  distance  from  the  planet’s 
surface.  Solutions  are  developed  for  the  non-rotating  Earth 
equations  of  motion  and  for  one  of  the  rotating  Earth 
equations  of  motion  using  the  method  of  matched  asymptotic 
expansions . 

It  is  shown  that  the  non-rotating  Earth  assumption 
produces  incorrect  entry  trajectory  results  for  entry 
orbital  inclination  angles  between  0.5  and  75.0  degrees  and 
vehicle  speeds  ranging  from  circular  orbital  velocities  to 
low  supersonic  speeds.  However,  a  variety  of  realistic 
trajectory  states  exist  where  some  of  the  non-rotating  Earth 


equations  of  motion  are  valid.  Three  of  the  non-rotating 
equations  of  motion  are  found  to  be  valid  for  the  same  entry 
trajectory  states.  Other,  independent  trajectory  states 
exist  where  a  fourth  non-rotating  Earth  equation  of  motion 
is  valid.  A  fifth  equation  of  motion  is  never  valid  for  the 
ranges  of  orbital  inclination  angle  and  speeds  investigated. 
Trends  in  the  results  of  the  trajectory  states  of  validity 
are  discussed  and  methods  to  estimate  some  of  these  states 
are  presented . 


INVESTIGATION  OF  THE  VALIDITY  OF  THE 
NON-ROTATING  PLANET  ASSUMPTION  FOR 
THREE-DIMENSIONAL  EARTH  ATMOSPHERIC  ENTRY 


I.  Introduction 


Analytical  studies  often  have  the  advantage  of 
displaying  valuable  solution  trends,  giving  more  insight  to 
the  problem  and  its  solutions  than  corresponding  numerical 
work.  Simple  and  accurate  analytical  methods  to  find 
solutions  to  the  equations  of  motion  for  high  speed  flight 
vehicles  supplement  more  complex  and  unwieldy  numerical 
methods.  In  past  analytical  work  on  lifting  atmospheric 
entry,  the  limiting  assumptions  of  planar  entry  and  a  non- 
rotating  planet  were  common.  The  equations  of  motion  for 
planar  atmospheric  entry  of  a  non-rotating  planet  are  highly 
nonlinear;  adding  rotating  planet  and  non-planar  effects  to 
the  equations  of  motion  make  them  even  more  difficult  to 
solve  analytically.  Hence,  the  current  literature  has  no 
investigations  which  analytically  solves  the  equations  of 
motion  for  three-dimensional,  lifting  atmospheric  entry  of  a 
rotating  planet. 

The  Non-Rotatine  Planet  Assumption 


Although  common  in  most  analytical  entry  trajectory 
analyses,  the  assumption  of  a  non-rotating  planet  model  can 
produce  significant  errors  in  trajectory  results.  Since 


most  planets  have  a  significant  rotation  rate,  the  rotating 
planet  assumption  will  give  more  accurate  entry  vehicle 
performance  results,  especially  for  lifting  vehicle  range 
and  time  of  flight  calculations.  Inherent  in  the  concept  of 
a  rotating  planet  for  atmospheric  entry  is  the  assumption 
that  the  planet’s  atmosphere  rotates  with  the  planet  at  a 
constant  rate.  This  assumption  is  more  accurate  than  the 
non-rotating  planet/atmosphere  assumption  but  is  still  not 
ideal .  The  atmosphere  of  a  rotating  planet  can  be  viewed  as 
a  boundary  layer  with  rotation  rates  which  change  with 
altitude.  Near  the  planet’s  surface  the  atmosphere  rotates 
at  about  the  same  rate  as  the  planet.  As  altitude 
increases,  the  atmosphere  rotates  with  a  decreasing  rate, 
and  at  high  altitudes  this  rate  is  near  zero.  Hence,  the 
true  effect  of  a  rotating  atmosphere  is  therefore  bounded  on 
one  end  by  the  non-rotating  planet  solutions  and  on  the 
other  end  by  the  rotating  planet  solutions.  In  this  study 
is  is  assumed  that  the  rotating  planet  solutions  are  ideally 
valid . 

The  maximum  effect  of  the  rotating  atmosphere  (Chapman, 
1959:5)  on  a  flight  vehicle  is  easily  calculated  for  a 
spherical  planet.  This  maximum  occurs  at  the  equator  for 
minimum  altitude  and  is  given  by  the  ratio  of  the  planet’s 
angular  velocity  to  the  circular  orbital  velocity  at  the 
planet’s  surface.  The  planet’s  rotational  velocity  on  the 
equator  is  given  by 


where  u  is  the  planet’s  rotation  rate  and  r  is  the  radius 

of  the  planet.  The  circular  orbital  velocity  at  the  surface 
of  the  planet  is 

V  =  (g  r  )5  /2 

where  g^  is  the  gravitational  acceleration  at  the  surface  of 
the  planet.  The  ratio  is  given  by 


This  ratio  gives  an  indication  of  the  possible  error 
introduced  to  velocity  calculations  for  a  non-rotating 
planet  model.  For  Earth,  this  ratio  is  approximately  0.06  . 
For  Mars,  Jupiter,  Saturn,  and  Venus,  this  ratio  is 
approximately  0.07,  0.30,  0.40,  and  0.0,  respectively  (Vinh 
and  others,  1980:3).  Hence,  for  near-equatorial  atmospheric 
entry,  the  maximum  velocity  error  is  negligible  for  Venus, 
significant  for  Earth  and  Mars,  and  very  large  for  Jupiter 
and  Saturn.  The  velocity  error  associated  with  the  non¬ 
rotating  planet  assumption  causes  an  even  larger  error  in 
the  calculation  of  aerodynamic  forces.  Since  lift  and  drag 
are  proportional  to  the  square  of  velocity  (as  seen  in 
Section  II),  the  error  introduced  to  the  aerodynamic  force 
calculations  can  be  a  maximum  of  about  0.14,  0.12,  0.60,  and 
0.80,  for  Mars,  Earth,  Jupiter,  and  Saturn,  respectively. 


A  recent  study  of  AOTV  (Aero-assi sted  Orbital  Transfer 
Vehicles)  trajectories  about  Earth  investigated  possible 
trajectory  simulation  error  due  to  the  non-rotating  planet 
assumption.  In  this  study  it  was  found  that  the  non- 
rotating  planet  model  caused  velocity  errors  which  gave 
dynamic  pressure  differences  ranging  up  to  10  -  14% 

These  differences  may  cause  underprediction  of  the  final 
altitude  and  overprediction  of  the  attainable  orbital 
inclination  change  in  a  non-rotating  Earth  analysis.  It  was 
concluded  that  rotating  Earth  effects  must  be  included  for 
realistic  AOTV  trajectory  simulation  (Ikawa,  1986:1,9). 

Another  recent  study  (Miller,  1986:14)  noted  significant 
errors  in  values  of  range  and  time  of  flight  for  equilibrium 
glide  entry  trajectories  when  the  Earth’s  rotation  was 
neglected.  Differences  in  trajectory  results  between  the 
rotating  and  non-rotating  cases  were  found  to  be 
significant,  especially  for  trajectories  beginning  at  speeds 
greater  than  15,000  ft/s. 

Vinh  (Vinh  and  others,  1980:3)  argues  that  inclusion  of 
a  rotating  atmosphere  into  an  analytical  entry  study  would 
cause  unwarranted  overcomplication  to  the  problem,  possibly 
overshadowing  the  advantages  of  an  analytical  versus 
numerical  analysis.  However,  for  many  atmospheric  entry 
trajectories,  such  as  multiple  pass  aerobraking,  synergistic 
orbital  plane  change,  and  general  high  L/D  vehicle 
trajectories  with  shallow  entry  flight  paths,  trends  in  the 
solutions  caused  by  the  rotating  planet  and  its  atmosphere 


may  be  important.  For  these  types  of  trajectories,  the 
error  in  the  calculation  of  the  aerodynamic  forces  is  more 
prominent  due  to  the  relatively  large  flight  times  within 
the  sensible  atmosphere. 


Problem 

Because  of  their  complexity,  the  equations  of  motion  for 
three-dimensional,  lifting  atmospheric  entry  of  a  rotating 
planet  have  not  been  analytically  solved.  However,  the 
assumption  of  a  non-rotating  planet,  common  in  most 
analytical  entry  analyses,  has  been  shown  to  produce 

i 

significant  errors  in  some  solutions  for  the  lifting  ( 

atmospheric  entry  of  Earth.  An  investigation  of  the  general 

validity  of  the  non-rotating  planet  assumption  for  three- 

dimensional  Earth  atmospheric  entry  is  needed.  In  addition, 

the  existence  of  trajectory  states  where  the  rotating  planet 

terms  in  the  equations  of  motion  are  negligible  should  be 

investigated.  This  would  indicate  the  existence  of 

trajectory  states  where  existing  solutions  to  the  non- 

rotating  equations  of  motion  are  valid  for  Earth  lifting 

atmospheric  entry. 

Scope 

In  this  effort,  the  three-dimensional  exact  equations  of 
motion  for  lifting  atmospheric  entry  are  expanded  to  include 
a  rotating  planet  model .  The  rotating  planet  terms  in  the 
equations  of  motion  for  Earth  atmospheric  entry  are 
examined.  Solutions  are  developed  for  the  non-rotating 


equations  of  motion  and  for  one  of  the  rotating  Earth 
equations  of  motion.  This  was  accomplished  by  treating 
atmospheric  entry  as  a  boundary  layer  problem,  and  applying 
method  of  directly  matched  asymptotic  expansions.  A  variety 
of  realistic  Earth  entry  trajectory  states  are  shown  to 
exist  where  some  of  the  non-rotating  equations  of  motion  are 
valid  for  a  rotating  Earth.  This  validity  is  coordinate 
dependent  since  sigularities  exist  in  the  equations  of 
motion.  Entry  trajectory  state  examination  is  limited  to 
orbital  inclination  angles  between  0.5  and  75.0  degrees, 
where  most  Earth  atmospheric  entry  occurs,  and  vehicle 
speeds  ranging  from  circular  orbital  velocity  to  low 
supersonic  speeds  where  terminal  maneuvers,  such  as  landing 
approaches,  are  usually  initiated. 


Assumptions 

The  planet  is  modelled  in  this  analysis  by  a  sphere 
having  a  central  gravitational  force  field  obeying  the 
inverse  square  law.  A  strictly  exponential  atmosphere, 
rotating  at  the  same  rate  as  the  planet,  is  assumed  with 
density  as  a  function  of  radial  distance  from  the  planet’s 
surface.  In  this  effort,  the  only  forces  acting  on  the 
lifting  vehicle  are  assumed  to  be  gravity,  lift,  and  drag; 
magnetic,  solar  wind,  and  other  forces  are  considered 
negligible.  The  lifting  entry  vehicle  is  modelled  as  a 
point  mass  in  a  three  degrees-of-f reedom  analysis.  Constant 
lift-to-drag  ratio  is  assumed  along  the  trajectory  and  a 


ballistic  coefficient  is  specified  for  each  flight 
vehicle/atmosphere  under  study.  Angle  of  attack  was  not 
explicitly  modelled.  More  detailed  discussion  of  the 
approximations  and  assumptions  will  be  presented  in  Section 
III. 

Approach 

In  Section  II,  the  equations  of  motion  for  three- 
dimensional  lifting  entry  for  a  spherical,  rotating  planet 
are  derived.  In  Section  III  these  equations  of  motion  are 
transformed  into  a  form  more  convenient  to  examine  and 
solve.  The  equations  are  also  made  non-dimensional  and  a 
coordinate  transformation  is  undertaken.  In  Section  IV  the 
rotating  planet  terms  in  each  of  the  five  equations  of 
motion  are  examined.  These  terms  then  are  set  equal  to  zero 
and  checked  for  the  existence  of  real  solutions.  It  is 
shown  that  three  of  the  equations  of  motion  have  identical 
solutions  for  these  rotating  terms.  It  is  also  shown  that 
real  solutions  do  not  exist  for  the  rotating  terms  in  one  of 
the  equations  of  motion.  The  solution  to  the  equation  of 
motion  containing  these  rotating  terms  is  developed  in 
Section  V  along  with  the  solutions  to  the  non-rotating  Earth 
equations  of  motion.  These  solutions  are  derived  from  the 
rotating  Earth  equations  of  motion  using  the  method  of 
matched  asymptotic  expansions.  In  Section  VI,  the  solutions 


to  the  rotating  term  equations  in  the  other  four  equations 
of  motion  are  examined  in  more  detail.  Trajectory  states 


are  presented  where  some  of  the  non-rotating  Earth  equations 
of  motion  are  independently  valid  for  a  rotating  Earth. 
Methods  are  given  to  estimate  solutions  where  four  of  the 
non-rotating  equations  of  motion  are  valid  for  rotating 
Earth  entry.  Conclusions  and  recommendations  for  further 
study  are  presented  in  Section  VII. 


In  this  section,  the  equations  of  motion  are  derived  for 
three-dimensional  lifting  entry  of  a  rotating  planet.  A 
spherical,  rotating  planet  model  is  employed  and  it  is 
assumed  that  the  atmosphere  rotates  at  the  same  rate  as  the 
planet  with  rotation  rate,  u .  The  lifting  entry  vehicle  is 
modelled  as  a  point  mass  in  a  three  degrees-of-f reedom 
analysis.  Gravity,  lift,  and  drag  are  assumed  to  be  the 
only  forces  acting  on  the  vehicle;  magnetic,  solar  wind,  and 
other  forces  are  assumed  to  be  negligable.  Further 
discussion  of  assumptions  and  approximations  is  presented  in 
Section  III. 

Coordinate  Systems 

Figure  1  defines  the  planet  centered  coordinate  systems 
used  in  this  analysis.  The  planet’s  inertial  reference 
frame  has  coordinates  X,  Y,  and  Z,  with  unit  vectors  I,  J, 
and  K,  respectively.  The  X  and  Y  axes  lie  in  the  planet’s 
equatorial  plane  and  the  planet  rotates  about  the  Z  axis. 

The  rotating  planet-fixed  coordinate  system  has  axes  XQ ,  YQ , 

and  ZQ ,  with  unit  vectors  IQ ,  JQ ,  and  Kc ,  respectively. 
Another  view  of  the  planet  centered  coordinate  systems  is 
illustrated  in  Figure  2.  a  is  a  unit  vector  in  the  XY 


plane,  pointing  radially  away  from  the  planet.  6  is  defined 


Figure  2.  Coordinate  Systems  Revisited 


as  the  longitude  and  is  measured  from  the  XQ  axis  which 

rotates  about  the  KQ  direction  with  rate  u.  f  is  latitude 
and  is  measured  positive  from  the  equator  to  the  pole  in  the 
Kq  direction.  7  is  the  flight  path  angle  and  is  measured 

"positive  up"  from  the  local  horizontal  to  the  velocity 
vector.  The  heading  angle,  is  measured  from  a  constant 
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systems.  The  velocity  of  the  vehicle  with  respect  to  the 
inertial  frame,  X,  Y,  Z,  is  defined  as  the  sum  of  the 
velocity  of  the  vehicle  with  respect  to  the  rotating  frame 
and  the  cross  product  of  the  angular  velocity  of  the 
rotating  frame  and  the  radius  vector. 


IJK 


dr 

dr 

dt 

IJK  dt 

+  fi  x  r 

lift 

(2.10 


where  r  is  the  radius  vector,  extending  from  the  planet’s 

center  to  the  flight  vehicle,  ^  is  the  velocity  vector,  and 

^  is  the  angular  velocity  of  the  rotating  frame. 

The  inertial  acceleration,  a,  is  defined  as  the 
derivative  of  the  inertial  velocity. 


a 

d2? 

d2? 

+  ft  x  ^ 

IJK  dt2 

IJK 

dt2 

♦  *  f  dt/ 

ijk 

+  (ft  x  ft)  x  r 


dt  dr 

+  n  x  -T- 

*  f  U 

ljk 


7  1  < 

1J 


dfi  -► 
+  x  r 
C  if  Ub 

ijk 

+  i)  x  in  x  r 


ij 


This  can  be  simplified  to  the  following 


where 


IJK 


is  the  acceleration  of  the  vehicle  with  respect 


to  the  inertial  planet  centered  coordinate 
system . 


is  the  acceleration  of  the  vehicle  with  respect 


to  the  radius  vector, 
is  the  angular  velocity 

a  point  mass  vehicle  in 
planet,  is  equal  to  a 


of  the  rotating  frame.  For 

flight  over  a  rotating 
constant,  the  planet’s 


rotation  rate,  u 


is  the  Coriolis  acceleration. 


a  x 


is  the  Centripetal  or  Transport 
acceleration . 


To  apply  Eq  (2.12)  in  the  derivation  of  the  equations  of 
motion,  transformations  between  the  various  coordinate 
systems  are  required.  The  coordinate  transformation  from 
the  inertial  system,  XYZ,  to  the  vehicle  centered  system, 
xyz,  is  not  difficult  but  is  prone  to  algebraic  error 
because  of  the  many  intermediate  steps  and  variables 
involved.  For  ease  in  derivation,  this  transformation  of 
coordinate  systems  is  divided  into  a  number  of  simple  angle 
rotations . 
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1st  Rotation,  Inertial  to  Rotating.  The  first  rotation 
is  the  planet  rotation  about  the  Z  axis.  The  rotation  angle 
at  time  t  is  wt.  From  Figure  1  it  can  be  seen  that  the 
following  expressions  apply: 


IQ  =  Icoswt  +  Jsinwt 
J  =  -Isinwt  +  Jcoswt 


K0  =  K 

In  matrix  form 


lo 

J 

coswt  sinwt  0 

-sinwt  coswt  0 

I 

J 

.  Ko. 

0  0  1 

K 

and 


I 

coswt  -sinwt  0 

lo 

J 

= 

sinwt  coswt  0 

lo 

K 

0  0  1 

Ko 

(2.1) 


(2.2) 


2nd  Rotation,  Longitude.  The  second  rotation  (Figure  5) 
is  a  longitude  angle  rotation  about  the  Z  axis: 


a  =  I  cos©  +  J  sin©  and  j  =  -I^sin0  +  J^cos0 

o  o  o  o 


3rd  Rotation,  Latitude.  The  third  rotation  is  a 
latitude  angle  rotation  about  the  y  axis.  From  Figure  6 


I  =  acos^  +  K  sin^ 


and 


Jc  =  -asin^  +  K  cos 0 


The  results  from  the  longitude  and  latitude  rotations  are 


vv 


combined  and  presented  below: 


J 

k 


and 


cos^cos#  cos^sin #  sirup 

?o 

-sin#  cos#  0 

-cos#sin0  -sin#sin^  cos <p 

Ko 

cos<pcos6 

-sin# 

sin^cos# 

♦ 

1 

-sin6cos<p 

cos# 

sin^sin# 

* 

3 

-sirup 

0 

cos^ 

k 

(2.3) 


(2.4) 


At  this  point  enough  information  is  known  to  make 

coordinate  transformations  between  the  XYZ  inertial  frame 

and  the  rotating  xyz  vehicle  point  mass  frame.  However,  the 

transformation  between  the  xyz  and  the  X0Y0Z0  frames  is 

required.  Figure  4  presented  the  vehicle  centered 

coordinate  systems  that  are  used  in  the  equations  of  motion. 

6  is  a  unit  vector  that  is  used  to  make  the  intermediate 

coordinate  transformations  easier  to  follow. 

4th  Rotation,  Heading  Angle.  The  fourth  rotation  is  a 
heading  angle  rotation  about  the  x  axis.  From  Figure  7 


B  =  jcosj}  +  ksin^fr  and  k0  =  -jsin^  +  Jccos^ 

5th  Rotation,  Flight  Path  Angle.  The  fifth  rotation  is 
a  flight  path  angle  rotation  about  the  zQ  axis.  It  can  be 

seen  from  Figure  8  that 
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loos')  -  Bsin^ 


and 


iBin')  +  £>cos7 


xo  = 


3n= 


The  results  from  the  heading  angle  and  flight  path  angle 
rotations  are  combined  and  presented  below: 


and 


COS') 

sin') 

0 


-sin^cos^  -sin^sin^ 
cos'ycos  y  cos^sin  y 

-siny  cos  y 


cos^cob  6  -sin') 

sin^cos^i  cos^cos  y 
sin^sinif)  -cos^sin  y 


0 

simp 
cos  if) 


3 


o 

Jo 

1L 


(2  .  b 


(2.6 


Combining  all  five  rotations  gives  the  relations  for  a 
complete  coordinate  transformation  between  the  x0y0z0  and 

XYZ  systems . 


COS') 

-sin^cosy* 

-sin7siny 

cosycos©  cosysin0 

siny 

-» 

Jo 

= 

sin') 

cos^cosy 

cos7sin^ 

-sin0  cos0 

0 

kJ 

0 

-siny 

cosy 

-cosfisiny  -sin^siny 

cosy 

r 

i  r  ~ 

1 

coswt  sinwt  0 
-sinu/t  coswt  0 


I 

J 


(2. 


where  c  is  a  unit  vector  in  the  direction  of  the  lift  force. 
From  Figure  6  it  can  be  seen  that  aerodynamic,  lift  acts 

in  the  c  direction,  drag  acts  in  the  -j  direction  and  the 

vehicle’s  velocity  vector  points  in  the  +j  direction. 


Hence,  lift  is  given  by 

£<  =  Lc  =  (Lcosa)I0+  (Lsina)k0 

From  Eqs  (2.9)  and  (2.5) 


£  =  Lcoscr  ^cos'll  -  sin'ycos^j  -  sin^sin^Jcj 
+  Lsincr  [-sin^j  +  cos  *t] 


(2.13) 


£.  =  L (cosacos7)  J  -  L (cosasin^cos^  +  sinasin\t)j 

+  L(sinacos^  -  cosasin^sin^) lc  (2.14) 

The  drag  components  are  found  in  terms  of  the  vehicle 
centered  unit  vectors  by  application  of  Eq  (2.6): 

S  =  -Djq  =  -D(sin7)i  -  D(cos7cos^)j  -  D (cos7sin^) Jc  (2.15) 

The  vehicle’s  local  or  relative  velocity  with  respect  to  its 
own  reference  frame  is 


V(sin7)i  +  V(cos7cos^)j  +  V(cos7sin^)  lc 


(2.16) 


The  planet’s  rotation  needs  to  be  accounted  for  to 
obtain  the  vehicle’s  inertial  velocity  from  Eqn  (2.16).  A 
velocity  component  due  to  the  planet’s  rotation  is  added  to 


the  vehicle’s  local  velocity  to  form  the  inertial  vehicle 
velocity.  The  planet’s  rotation  velocity  component  is  a 
function  of  latitude;  at  the  equator  this  velocity  component 
has  its  maximum  value  and  at  the  poles  it  is  zero. 

Vehicle’s  Inertial  Velocity  =  P  - 

IJK 

Velocity  component  due  to  planet  rotation  =  ^ 

Vehicle’s  local  velocity  = 

ijt 

$  ~~~  =  $  +  $  (2.17) 

IJK  ijk  r 

=  Vrj  =  wrcos^j  (2.18) 

Combining  Eqs  (2.16),  (2.17),  and  (2.18)  gives  an  expression 

for  the  inertial  velocity  of  the  flight  vehicle. 

^ - =  V(sin7)i  +  V(cos^cos^)  +  wrcos^)  j 

IJK 

+  V(cos7sin^i)  k  (2.19) 

Another  expression  for  inertial  velocity  can  be  derived 
from  Eq  (2.10)  and  compared  to  Eq  (2.19)  to  produce  three  of 
the  equations  of  motion. 

A  dr  dr  *  „  -» 

IJK  dt  IJK  dt  Ijfc  I3R 

In  this  equation  tl  is  the  vector  sum  of  two  rotation 

rates  between  three  coordinate  systems,  the  XYZ ,  X  Y  Z  ,  and 

J  ’  0  0  0 


xyz  systems.  Therefore,  fl  is  the  sum  of  the  angular 
velocity  of  the  rotating  planet  frame  about  the  inertial 
frame  and  the  angular  velocity  of  the  vehicle  centered 
reference  frame  about  the  rotating  planet  frame.  This  can 
be  expressed  by 


a  =  a 


XYZ— »X  Y  Z 

0  0  0 


XYZ 

0  0  0 


►xyz 


From  Figure  1 


-  £  +  i  +  $  -  [w  +  H]k0  -  If  5 


and  with  Eqn  (2.4) 


=  [w  +  <it]sin^  -  It3  +  [u  +  ff]cos^ 


(2.20) 


Since 


ri 


dr 

dt 


m 


dr?  dl  _  dr? 

dt1  +  rdt  »•*«■  dt1 


and 


T[w  +  af]cos^3  +  It6 


(2.21) 


Hence,  the  second  equation  for  inertial  velocity  in  terms  of 
the  vehicle  centered  reference  frame  coordinates  is 


=  If1  +  rcos^(w  +  |£)j  +  r|£fc 


(2.22) 


Three  kinematic  equations  of  motion  are  derived  by 
equating  the  two  inertial  velocity  expressions,  Eqs  (2.19) 
and  (2.22),  and  comparing  like  terms. 


+  rcos <p(w  +  4£)  j  +  =  V(sin7)i  + 


+  V(cos7cos rf>  +  urcos^)  ]  +  V(cos7sin^)  k 
Kinematic  Equations  of  Motion. 


i  terms : 


St  =  Vsin^ 


(2.23 


j  terms:  -r- 


VcosTfcos; 

rcos^ 


(2.24 


k  tei 


Vcosnfsin! 


(2.25 


Derivation  of  the  Force  Equations  of  Motion 

To  derive  additional  equations  of  motion,  the  inertial 
acceleration  is  calculated.  From  Eq  (2.12) 


=  a  +  2fl  x  r 

ilk 


+  x  r  +  ft  x  fixr 

ijk  ljk  ljk 


Each  of  the  four  terms  in  the  above  equation  is  derived 
separably  for  clarity. 


1st  Acceleration  Term.  a 


From  Eqs  (2.21)  and  (2.23) 


Therefore  a  =  si 

ijfc  L 


in7^  +  Vcos75?J  i 


2nd  Acceleration  Term.  X  r 


From  Eqs  (2.20) 


j*  d01  .  ,*  d  <t>->  d0]  ,r 

n  =  u  +  smpi  -  +  u  +  cospk 


and  from  above 


-*  dr » 

r  *•»(•  _  dt1 
ijk 


dr?  ,r  .  ? 

=  TT1  =  Vsin7i 


Taking  twice  the  value  of  the  cross  product  of  these  two 
equations  gives  the  second  inertial  acceleration  term. 

2?l  x  r  =  |2Vsin7cos^  (w  +  ^)]j  +  2Vsin7^jk  (2 

ljk 


3rd  Acceleration  Term.  ^  x  r 


Taking  the  time  derivative  of  Eq  (2.20)  gives 


*  r  del  ,d^r  d7e  .  ,?  .  r  .  dei  .  ,di  d7<t> , 

=  lw +  dtj^dt1  +  +  r  +  dtjsin^dt  - 


■  [w  +  Jt]sin4tfc  +  TT006^  +  ["  +  5t]cos^H 

at 


Rewritten 


d  e 


Therefore,  the  third  term  of  Eq  (2.12) 


4th  Acceleration  Term.  fl  X  fl  X  r 


Taking  the  cross  product  of  Eqs  (2.20)  and  (2.21)  gives 


X  ?  =  r[w  +  ff]cos^  +  rft* 

1  J  K 


Taking  another  cross  product  produces  the  fourth  inertial 


acceleration  term. 


ft  x  (ft  x  r  ) 

*56 


r  d01 2  2  . 

r[w  +  dtj  cos  * 


d2  ^  * 
+  r — l 


rft[W  +  ff]  sin^  3  +  r[u>  +  sin^cos^  k 


Total  Inertial  Acceleration.  Combining  Eqs  (2.26), 


(2.27),  (2.29),  and  (2.30)  gives  the  total  inertial 


acceleration  in  terms  of  the  vehicle  centered  direction  unit 


vectors . 


a - =  Is 

IJK  l 


.  dv  v  di  r  Aey  2 .  a 2  4  * 
in7dt  +  Vcos7dt  -  rr  +  dtj  cos  *  -  r7T  1 

at/  . 


SKM 


+  2Vsin7cos^w  +  ^-J  +  r — —cos ^  -  2]u)  +  j^-Jrsi 

L  dt 


r[w  +  sin^cos^  +  r— ^  +  2Vsin7^  k 


(2.31) 


This  equation  can  be  simplified  by  substitution  of  the 
three  kinematic  equations  of  motion  that  were  previously 
derived,  Eqs  (2.23),  (2.24),  and  (2.25).  Because  of  the 
size  of  this  new  equation,  each  direction  component  of  the 
inertial  acceleration,  Eq  (2.31),  is  examined  separably 
below . 

Ith  Inertial  Acceleration  Component.  The  Ith  term  of  Eq 
(2.31)  is 

.  dV  d7  T  del2  2  .  d 7<b 

ith  =  sm7^  +  Vcos7j£  -  r[w  +  cos  0  -  r~f 

at  . 

Substituting  in  for  ^  and  ^  with  Eqs  (2.24)  and  (2.25) 
Ith  =  sin7j£  +  Vcos7jJ  -  rw2  cos2  0  -  2u;r  [Vc°^°s^]  cos 2  <p 


o  1/2  2  2 

2  .  V  cos  nrcos 


,,2  2.2 
V  cos  nrsin 


r  cos  f 
This  can  be  simplified  to 


dV  d'l  V2  2 

=  sin7^  +  Vcos7^f  -  —'cos  ^  -  2u/Vcos7cos^cos0 


2  2 

u  rcos  <p 


This  equation  simplifies  to  the  following  expression 


2 

j  th  -  ^C°cosd)S^  TCOS(P  +  2wVsin7Cos <p  +  2^’ sin7cos^cos if) 


2 

-  2wVsin^cos7sin|t  -  2—  ’tan^cos  7sin^cos ^ 


The  first  term  of  this  equation  is 


d  Ycoe'ycosit  ,  ,dV  .  ,d7 

dt  "  rcos 4>  TCOS<f>  =  cos7cosV'dt  "  Vsin7cos^ 


,dtf  VcoeTfcostf  ,dr  .  ,d  <t> 

-  Vcos7sin^ - r  2/  cos*dt  "  rsin*dt 

r  cos  q> 


Simplifying  and  substituting  in  Eqs  (2.23)  and  (2.25)  gives 
3t  VCrcosiS*  rcOS<f>  =  cos7cos ^  -  Vsin7cos^  -  Vcos7sin^ 


V2  V*  2 

-  -■  sin^cos'ycos^  +  -'cos  ^cos^sinjttan^ 


Therefore,  the  entire  j th  term  of  Eq  (2.31)  is 
jth  =  cos7cos^  -  VsinTfcos^^  -  VcosTfsin^g^ 

+  2wV^sin'Ycos0  -  sin^cos7sin^J 


s^[si 


—  *  cos'jcosf  |^sin7  -  COS7SI 


in^tan^j 


(2.33) 


kth  Inertial  Acceleration  Component.  The  kth  term  of  Eq 


(2.31)  is 


lcth  =  r^  +  2Vsin7^^  +  r  sin^cos^ 

The  first  term  in  this  fcth  component  requires  some 
manipulation.  Differentiating  Eq  (2.25)  and  multiplying  by 
r  gives 


d2d>  d 
r — =  r  •  — 


fcosTfsinii  .  ,dV  ,,  .  .  ,d7 

^ - *■  =  cos^smy^-  -  Vsin^siny^^ 


+  Vcos^cos^^  -  sin7cos7sin^ 


(2.34) 


Substituting  Eqs  (2.34),  (2.23),  and  (2.24)  into  the  lcth 
term  produces 


,V2 


2  2  o 

V  cos  nrcos  ifrsinij 


f-  .  1  _V.  .,  I  .  ,  ,  V  COS  T  cos 

kth  =  2—  •  sin'jcos'jsxnf  +  w  rsinpcosp  +  - rcosft — 

dV  d'Y 

+  2wVcos7cos^sinrf  +  cos'isinit-rr  -  ’Vein'isiurb-rr 


+  Vcos'ycos^Ji  _  Bin'jcos'jsinip 


This  equation  can  be  simplified  to  the  following  expression 


\r  2  dV 

kth  =  -■  COS7 (COS7COS  jttan^  +  ein^sini/i)  +  cos^sin^j^ 

-  Vsin^sin^^  +  Vcos'ycos^^  +  2u/Vcos7cos^sin0 


+  u 1  rsin^cos^ 


(2.35 


Total  Inertial  Acceleration.  Adding  together  the 
ith,  j  th ,  and  kth  acceleration  terms  given  by  Eqs  (2.32), 
(2.33),  and  (2.35)  gives  an  expression  for  the  total 
inertial  acceleration  on  the  vehicle. 

-»  dV  d'Y  V2  2 

a - =  sin'll  +  Vcos'f^f-  -  2uVcos'7cos^cos^>  -  —  *  cos  7 

-  w2  rcos ^  <p  1  +  cos7cos^^-  -  Vsin7cos^^  -  Vcos^sin^^^ 
+  2Vw  (cos<psinr)  -  sin^cos^sin ifi) 


+  — • cos^cos^ (sin^  -  sin^cos^tan^)  3 


•  j.dV  v  .dtf 

+  cos^sin^  +  Vcos^cosy^-^ 


-  Vsin7siny/^  +  2u/Vsin^cos7cos^ 


(2.36 


2  *  V  .  2  c 

+  u  rsin^cos^  +  — * C0S7 (sin7sin^  +  cos  ^tan^cos7)  k 


Another  equation  for  inertial  acceleration  can  be 
derived  by  examining  the  forces  on  the  vehicle.  It  is 
assumed  that  the  only  forces  acting  on  the  vehicle  are 


gravity  and  aerodynamic  lift  and  drag. 


Therefore 


and 


IJK 


d$ 

dt 


IJK 


=  ;(£  -  6) 


+  e 


(2.37) 


(2.38) 


where  the  gravitational  force,  g,  is  a  function  of  r  and 
acts  in  the  negative  radial  direction. 


g  =  ~g(r)i  (2.39) 

Lift  and  drag  are  given  by  Eqs  (2.14)  and  (2.15): 

t.  =  (Lcosctcos7)  1  -  (Lcosasin7cos^  +  Lsinasin^)j 
-  (Lcosasin7sin^  -  Lsinacos^)lc 

5  =  (-Dsin7)i  -  (Dcos7cos^)  j  -  (Dcos7sin^)  1c 

With  substitution,  a  second  equation  is  found  for  the  total 
inertial  acceleration  of  the  vehicle: 


IJK 

[L  D  .  ]t 

— cosacos7  -  ^sm7  -  gl  1 

1 

[L  •  _  ^  .  L  . 

D 

-  ^cosasin7sin^  -  ^sinacosjj  +  j^cos7sin^J  &  (2.40) 


By  setting  the  two  different  equations  for  inertial 


acceleration  equal  and  comparing  like  terms,  three 
additional  equations  of  motion  can  be  derived.  Equating  Eqs 
(2.40)  and  (2.35)  produces  the  following  expressions: 


I  Direction  Term. 


dV 

sln7dt 


d'V  V  2 

Vcos^j-^  -  2u/Vcos7cosyicos^  -  —• cos  7 


2  2 

u  rcos  <p 


-(Lcosacos'j  -  Dsin7)  -  g 


(2.41) 


3  Direction  Term . 


cos 


dV  .  d7 

S7~p£  -  Vsinnf^f  - 


n7dt  "  Vcos7tan^at 


2Vw 


+  (cos^sin7~sin$cos7siny)  + 


V2 

—  *  COS7  (sin7-sin^cos7tan^>) 


-1 

m 


(Lcosasin7 


+  Lsinatan^  +  Dcos7)j 


(2.42) 


Jc  Direction  Term. 


dV 

cos7dt  + 


wCosj  dj£ 
tan^i  dt 


Vsin7^t  + 


2wVsi^c°s7 

tany 


+ 


w 


si 
r — 


n^cos^ 

sin^ 


V" 

+  —  •  COSTf 

r  ' 


sin7  + 


cos^tan^cosnf 

tan^i 


(Lcosasin7 


y  sing 
^tan^ 


+ 


DCOS7) 


(2.43) 


These  three  coupled  equations  can  be  reduced  by  some 
manipulation.  Multiplying  Eq  (2.43)  by  -1  and  adding  the 
product  to  Eq  (2.42)  gives 


-  Vcos7tan 


+ 


2Vw 

COBlfl 


(cos^sin7  - 


sin^cos7sin^) 


-  ?-cos7[sinycos7tan^]  -  VffjJJ  -  2^Vsig^°t 

_  w2  rsin^cos^  _  V2.cos  rcos^tan^cosj  ] 
sint t>  r  '  L  tanti  J 


cosytan^cos' 

ta,nip 


L  .  ,  ,  L  sina 

=  -  — sinatantt  -  —  •- - r 

m  T  m  tan rp 


By  combining  terms  and  noting  the  trigonometric 


identities 


2  2  _ 1  2  2 

1  +  tan  ^  =  (cos  ip)  and  sin  f  +  cos  y  =  1 

the  above  expression  can  be  rewritten  to  form  a  new  equation 
of  motion : 


L . sing 
m  cos7 


+  2Wui  (cos0tan7sin^  -  sin^) 


-  lj  r 


cos7 


-  — • cos7cos^tan^ 


(2.44) 


Substituting  Eq  (2.44)  back  into  Eq  (2.43)  produces 


dV  dT  V2 

0  =  cos7^  +  2Vu;cos0sin7cos^  -  Vsin7j^  +  -•  cos7sin7 

+  ^cosasin7  +  j^cos7  +  w2  rsin^cos^sin^ 


(2.45) 


Multiplying  Eq  (2.41)  by  -  and  adding  this 

product  to  Eq  (2.45)  gives  another  equation  of  motion: 


dtl51^  + 


cos2  Trl  V2  r cos2  nr  .  1 

sin7  J  r  '  lsin7  'J 


L  F  .  cos2 

+  — cosg  sin7  +  — : - 

m  l  '  sin7 


Ur1]  "  +  2Vu/cos^cos^[sin7 


Rewritten , 


?S  +  w2  rcosi  sin^sintf  +  cos^^?^] 
sin^  J  r  L  ~  ~  rsm7J 


=  ^*cos7  +  “cose  -  gcos^  +  2Vwcos0cos^ 
+  w2  rcos^  |sin^sin^sin,y  +  cos^cos^j 


(2.46) 


Substituting  Eq  (2.46)  into  Eq  (2.42)  produces  the  last 
equation  of  motion: 


dV  2  L  2  n\r  * 

sin7j^  +  —’cos  7  +  — c.os'jcosa  -  geos  7  +  2Vujcos')cos<pcosf 


u/2  rcos^cosip  cos0cos7  -  sin^sin^sin7j 


v  2  2  2  L 

-  2wVcos7cos^cos^  -  — *  cos  7  -  w  rcos  ^  -  -,cosacos7 


+  ^sin7  +  g 


This  can  be  rewritten 


dV  D  .  2  r  1 

=  -  —  -  gsin7  +  u  rcos0 I cos0sin7  -  sin^sin^cos7 1  (2.48) 


The  Equations  of  Motion 

In  summary,  the  following  six  equations  of  motion  for 
three-dimensional  atmospheric  entry  for  a  rotating  planet 
have  been  derived  and  are  listed  below  for  convenience. 

^  -  gsin7  +  w2rcos^fcos^sin7  -  sin^sin^cos7j  (2.49) 


=  —  'COS')  +  —  •  COSIT  -  geos'!  +  2Vwcos^cos^ 
+  w?  rcos^i  |sin^sin^sin')f  +  cos^cos'jfj 

vj*  =  ^  +  2Va/ |cos^tan7sin^  -  sin^>j 


u7  r 


sind>cos(ftcosV* 
cos' y 


-  -■  cos'ycos^tan^ 


dr 

dt 


-  Vsin7 


d6 

dt 


dt 


Vcos^cosip 
rcos  <p 

Vcos^sin  rj> 
r 


In  the  next  section,  these  equations  of  mot 


ion  are 


transformed  into  a  form  more  convenient  to  analyze. 


Ill .  Approximations  and  Manipulation 
of  the  Equations  of  Motion 


The  equations  of  motion  for  three-dimensional 
atmospheric  entry  for  a  rotating  spherical  planet  were 
derived  in  Section  II.  In  this  section,  assumptions  and 
approximations  used  in  this  thesis  are  defined  and 
discussed.  In  addition,  the  equations  of  motion  are 
transformed  into  a  form  more  convenient  to  examine  and  solve 
in  later  sections.  The  independent  variable  is  changed  from 
time  to  non-dimensional  altitude,  h,  and  the  equations  of 
motion  are  transformed  into  non-dimensional  form  by  the 
introduction  of  non-dimensional  variables.  A  coordinate 
system  transformation  is  undertaken  to  utilize  variables 
which  are  more  convenient  for  atmospheric  entry  analysis. 

The  equations  of  motion  derived  in  Section  II  are 
given  by 

dV  D  2 

-rr  =  -  —  -  gsin^  +  u rcos^ (cos^sin7  -  sin^sin^cos7)  (3.1) 


vdj  V2 
^dt  =  r  cos7 


+  u  rcos <p  ( 


ydj£  _  L.  sing 
dt  m  COS7 


^cosa  -  gcos7  +  2Vu/cos^cosyi 
sin^sin^sin7  +  cos^cos 7) 

V2 

-•  cos7cos|ftan^  +  2Vu/ (cos^tan7sin^ 


(3.2) 


sin^) 


0 

u  rsinicostfcosf 


dr  ,,  . 
dt  =  Vsin^ 


d 6  _  VcoBTfcos{ 
dt  rcosi 


d<t>  _  Vcos^fsin; 
dt  r 


(3.4) 


(3.5) 


(3.6) 


Note  that  the  equations  of  motion  for  three-dimensional, 
non-rotating  planetary  entry  can  be  derived  from  the 
rotating  equations  simply  by  setting  the  planet  rotation 
rate,  w,  to  zero.  The  equations  of  motion  for  the  non¬ 
rotating  planet  assumption  are  therefore 


dV  D 

dt  =  "  m  -  esm 7 


=  £ 


>7  +  -coscr  -  gcos7 


2 

,,dt£  L  sina  V  , 

y~rr  =  —  ■ - -  -  —  •  cosnrcosifltan d> 

dt  m  cosnr  r  i  r  r 


dr  . 

dt  =  Vsin? 

d0  _  Vcos7cos| 
dt  rcos<p 

d  d>  Vcos7sint 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


As  expected,  comparison  of  Eqs  (3.1)  -  (3.6)  and  Eqs 
(3.7)  -  (3.12)  shows  that  the  equations  of  motion  are 
significantly  more  complicated  when  the  earth’s  rotation  is 
accounted  for. 


A  few  more  variables  require  definition  at  this  point. 


Let  y  be  altitude  and  r+  be  average  planetary  equatorial 

radius.  Since  r  is  the  radius  measured  from  the  center  of 
the  spherical  planet  to  the  flight  vehicle,  h  is  defined  as 
the  non-dimensional  altitude  and  is  given  by 

h  =  (3-13) 

★ 


where 

r  =  r^+  y  =  r^+  hr^  =  r^(l+h)  (3.14) 

Therefore,  =  r+  (3.15) 


,  dt  _  dr/dh  _  r* 
dh  dr/dt  Vsin'y 

These  relations  will  be  used 
on  the  following  pages. 


(3.16) 


in  the  approximations  discussed 


Assumptions  and  Approximations 

The  rotating  planet  assumption  and  the  errors  associated 
with  the  non-rotating  approximation  were  discussed  in 
Section  I.  Further  assumptions  and  approximations  are 
presented  below. 

Spherical  Planet  Assumption.  The  approximation  of  an 
oblate  planet  by  a  sphere  is  very  common  in  analytical 
flight  mechanics  analyses  and  is  used  here.  It  is  a 
reasonable  assumption  for  planets  having  small  ellipticity 
such  as  Earth,  Venus,  and  Mars  but  may  not  be  as  reasonable 
for  planets  such  as  Jupiter  and  Saturn  which  have  relatively 


large  equatorial  bulges  and  ell iptici ties  a  magnitude 
greater  than  Earth’s  (Chapman,  1958:2).  For  either  case, 
error  introduced  by  the  spherical  planet  approximation  is 
generally  small  for  near  equatorial  entry  trajectories. 

Spherical  Atmosphere  Assumption.  The  spherical  planet 
approximation  leads  to  the  assumption  that  the  atmosphere  is 
spherically  symmetric  about  the  planet.  In  reality,  planets 
are  oblate  spheroids  causing  their  atmospheres  to  also  have 
an  oblate  form.  In  addition,  other  significant  deviations 
from  the  spherical  model  occur.  A  diurnal  density  bulge 
occurs  over  part  of  Earth’s  sunlit  side  due  to  solar 
heating.  Solar  storms  and  fluctuations  in  a  planet’s 
magnetic  field  can  cause  significant  changes  in  density  for 
a  given  altitude  (Wiesel,  1986:66-69).  However,  since  these 
effects  generally  occur  at  altitudes  where  aerodynamic 
forces  are  minimal  (and  inclusion  of  more  sophisticated 
density  models  may  introduce  overcomplication  of  the  entry 
problem)  these  effects  are  generally  assumed  to  be 
negligable.  The  approximation  of  a  spherical  planetary 
atmosphere  is  perhaps  the  most  limiting  assumption  to  the 
atmospheric  model.  This  approximation  is  better  for  the 
Terrestial  planets,  with  their  slow  rotation  rates,  than  for 
the  large,  outer  planets  (Duncan,  1962:276). 

Gravitational  Model .  The  spherical  planet  approximation 
also  leads  to  the  assumption  of  an  inverse  square 
gravitational  field.  This  is  given  by 


where  is  the  gravitational  acceleration  at  the  planet 
surface . 

Atmospheric  Density  Model .  As  previously  discussed,  the 
planet’s  atmosphere  is  assumed  to  be  spherical  and  to  rotate 
with  the  planet  with  a  constant  rotation  rate,  w .  The 
planet  atmosphere  is  mathematically  modelled  by  an 
exponential  atmosphere  with  the  inverse  atmospheric  scale 
height,  /J .  This  is  a  very  common  atmospheric  model  that  has 
been  successfully  utilized  in  many  studies  on  planetary 
atmospheric  entry.  Atmospheric  density,  p,  is  given  by 
(Chapman,  1959:4) 

^  =  -/Jdr  (3.18) 

P 

It  is  sometimes  assumed  that  the  product  of  the  inverse 
of  the  atmospheric  scale  height  and  the  vehicle’s  distance 
from  the  planet’s  center  is  constant  for  a  given  atmosphere. 
With  this  model 

^  r  or  p  =  (r/r^)  (3.19) 

where  p  is  the  density  at  the  surface  of  the  planet. 

The  product  fir  has  been  approximated  for  many  of  the 
planet  atmospheres.  It  has  values  on  the  order  of  1000  for 
most  of  the  planets;  the  mean  value  of  /Jr  is  approximately 


900  and  350  for  Earth  and  Mars,  respectively  (Vinh  and 
others,  1980:5). 

For  this  study  a  strictly  exponential  atmospheric  model 
is  employed  where  the  the  inverse  atmospheric  scale  height 
is  assumed  to  be  constant.  This  approximation  allows  the 
atmospheric  density  equation  to  be  written  in  the  form 

P  =  P*e~^y  =  P*e~h//e  (3.20 


where  e  is  a  small  number  given  by 


(3.21 


Note  that  is  dependent  on  what  planet  is  studied;  for 
Earth,  the  average  scale  height  is  about  7.1  kilometers 
(Vinh  and  others,  1980:5)  or  about  23,300  feet. 

For  planetary  atmospheric  entry,  r  is  approximately 
equal  to  r^ .  This  is  an  accurate  approximation  because  the 

thickness  of  an  atmosphere  is  generally  very  small  compared 
to  the  planet’s  radius.  For  example,  the  upper  altitude 
limit  of  Earth’s  sensible  atmosphere  is  often  taken  to  be 
350,000  feet.  This  value  is  only  about  1.795  of  the  Earth’s 
radius.  Approximating  r  by  r^  leads  to 


1  =  _1_ 

PT*  PT 

and  hence  the  values  given  in  the  literature  for  mean 
planetary  /Jr  are  considered  equivalent  to  f)r  .  e  is 


therefore  a  very  small  number  (approximately  equal  to  1/900 


for  Earth)  allowing  it  to  be  utilized  as  a  small  parameter 
in  the  asymptotic  expansions  of  Section  V.  Earth 
atmospheric  density  values  were  calculated  with  this  model 
and  plotted  in  Figure  10  with  values  obtained  from  the  1976 
U.S.  Standard  Atmosphere  (NOAA,  1976:Table  IV).  Comparison 
of  these  results  shows  that  the  accuracy  of  the  exponential 
density  model  is  reasonable  for  Earth  with  a  constant  value 
for  /Jr  of  900. 

Aerodynamic  Forces  and  the  Ballistic  Coefficient.  Lift 
and  drag  accelerations  are  given  by  the  following  familiar 
expressions : 


L  PSV*CL  ,  D  ^SV2CD 

m  “  2m  and  m  ~  2m 


(3.22) 


where  m  is  the  mass,  S  is  the  aerodynamic  reference  area,  V 


is  the  velocity,  C^  is  the  lift  coefficient,  and  C^  is  the 
coefficient  of  drag  for  the  flight  vehicle. 

The  non-dimensional  ballistic  coefficient,  B,  is  defined 
to  help  place  the  equations  of  motion  into  a  form  easier  to 
solve.  It  is  given  by 


P~SCT 


(3.23) 


and  is  specified  for  each  flight  vehicle  under 
consideration.  The  ballistic  coefficient  is  a  function  of 
the  vehicle’s  physical  characteristics  and  the  planet 
atmosphere  and  is  considered  to  be  a  constant  (Busemann  and 
others,  1976:18). 


'/>/ /  / 


cos'y  (1+h)  J 


Differentiating  Eq  (3.18)  with  respect  to  h  gives 


du  _  2V (1+h)  cos2  7 .  dV 


dh 


g  r 
°  *  * 


dh 


2  •  2  2 
2V  (1+h)  cosTfsinTr .  d^  V  cos  i 

gr  dh  g  r 

*  4c  4e 


(3.26) 


Substituting  Eqs  (3.20),  (3.23),  and  (3.25)  into  Eq  (3.22) 
gives  a  non-dimensional  equation  for  the  aerodynamic  drag 
acceleration : 


—  -  v2  — Qf  ^/e  -  fl rv5»  ^/e 


=  V  SCne 
m  2m  D 


=  pBV  e 


or 

D  ug^Bexp (-h/e ) 

m  e (1+h) cos2  7 

Likewise 

L  =  D.^L  =  CL  ug»Bexp (-h/e) 
m  m  Cjj  Cp  e  (1+h)  cos2  7 


(3.27) 


(3.28) 


Changing  the  Independent  Variable 

The  equations  of  motion  for  a  rotating  planet  are  now 
rewritten  using  the  expressions  derived  above.  To  change 
the  independent  variable  from  time,  t,  to  non-dimensional 
altitude,  h,  the  equations  of  motion  are  multiplied  by  Eq 
(3.16).  Direct  substitution  into  the  dV/dt  equation  gives 
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ug  r 


1  /  2 


• cos0cos7 (cos^sin^  -  sin0sin^cos7) 


is  equation  can  be  reduced  to 


Bexp (-h/e) 


ecos7sin7 
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COS7 
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tan7 


(l+h)3^3!1  /2 


u(l+h)J J 

(cos^isin7  -  sin0sin^cos7)  (3.29) 


Direct  substitution  into  the  d7/dt  equation  gives 


a_.. 


CL  ug»BexP (_h/e )  cosa  r^ (1+h) cos2 7 


cos 


€003731117 


(l+h)sin7  tan7  ug 


2  2 

u  r^(l+h)  sin^cos^cosf 

• (cos0tan7sin^-sin0)  -  - - - -  (3.31) 

u®* tan / 


The  d6/dt  equation  becomes 


dh  (1+h) cos^tan7 

and  the  dip  /  dt  equation  becomes 

dip  _  sin rp 

dh  (l+h)tan7 

The  dr/dt  equation  is  incorporated  into  the  other 


(3.32) 


(3.33) 


equations,  reducing  the  number  of  equations  of  motion  from 
six  to  five.  The  du/dh  equation  also  incorporates  the  dV/dh 
equation.  Combining  Eqns  (3.19)  and  (3.20)  allows  the  du/dh 
equation  to  be  rewritten  as 
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d 


2 


Substituting  in  for  the  dV/dh  and  d'y/dh  terms  using  Eqs 
(3.29)  and  (3.30)  eliminates  V  and  gives  the  du/dh  equation 
as  a  function  of  u. 
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This  rather  formidable  equation  reduces  to  the  following 
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The  du/dh  equation  can  be  further  reduced  to 
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(3.34) 


Another  variable  is  introduced  to  simplify  the  equations 
usemann  and  others,  1976:19). 


t  q  =  cos'y  h< 


dq  .  dnr 

dh  =  -Sin^dh 


(3.35) 


ing  these  relationships,  the  d'y/dh  equation  can  be 


written 


=  _ 3 _ [s2_  jl  _  ft. 

(l+h)lu  *J  CD 


cL  Be-h/€  rr,(l+h)l’/2 

~ - cosa  -  2qwcos$cosii  - 

CD  e  UZ* 


u?  r^q3  (1+h)  2  "^^(sin^sin^tan^  +  cos^) 


(3.36) 


The  equations  of  motion  for  three-dimensional,  rotating 
anetary  entry  have  now  been  transformed  from  Eqs  (3.1)  - 
.6)  to  the  following: 


2uBe 


-h/e 


(1+h)  esin7 


1  +  *tan7cosa  -  4w 


(1+h) ur 


• cosAcostf  -  2u 


2  <1+h>  r* 


>sin®smt 


(3.37) 


cosa  -  2qu/cos0cosd» 


(sin^sin^tan7  +  cos^) 


(3.38) 


-  w2  r  q3(l*h)2-£2^ 

V  /  Ug^ 


djt  _  ^L.Be  ^'//£sinq 
dh  -  Cp  ecos^sin^ 


cosTrcos^tan#  ^ 
(1+h) sin^ 


2u 

tan7 


(l  +  h)r  11/2 


u8. 


•  (cos^itan^sin^-sin^)  - 


2  2 

w  r^(l+h)  sin^cos^cos^ 
ug  tan7 


(3.39) 


dO  _  _ cos^ 

dh  (1+h) cos^tan7 


(3.40) 


dip  _  si  n^ 
dh  (l+h)tan7 


(3.41) 


The  Classical  Orbit.  Variables 

Up  to  this  point,  the  equations  of  motion  have  been 
presented  as  functions  of  latitude,  longitude,  heading 
angle,  and  other  variables.  This  form  has  been  often 
utilized  for  atmospheric  trajectory  simulation  by  numerical 
integration.  For  ease  in  studying  and  in  deriving  solutions 
in  later  sections,  these  equations  are  placed  in  terms  of 
the  classical  orbital  elements,  ft,  I,  and  a.  For  a  non¬ 
rotating,  spherical  planet,  these  variables  are  constants  of 
motion  for  non-atmospheric  flight.  This  characteristic 
greatly  simplifies  the  solution  derivation  for  the  non¬ 
rotating  planet  case.  I  is  defined  as  the  orbital 
inclination  angle,  fl  is  the  longitude  of  the  ascending  node, 
and  a  is  the  argument  of  latitude  at  epoch.  Basic  spherical 
trigonometric  relations  are  found  in  Appendix  A  and  these 
relations  are  applied  to  transform  the  variables  6,  $ ,  and  ^ 


51 


V.  V, 


in  terms  of  the  orbital  elements,  a,  fl,  and  I.  Figures  11 
and  12  show  the  geometry  of  the  two  sets  of  variables. 

The  following  relations,  derived  in  Appendix  A,  relate 
6 ,  <p ,  and  ip ,  and  a  ,  0 ,  and  I . 


sirup 

=  sinlsina 

(3.4 

cosl 

=  cosipcosip 

(3.4 

si  nip 

tan# 
tan  a 

(3.4 

sin  (0 

tanji 
”  tanl 

(3.4 

sin^ 

=  sinlcos(0-fi) 

(3.4 

COSQ 

=  cos#cos (0-fl) 

(3.4 

Differentiating  Eq  (3.43)  gives: 
sinl'dl  =  cos^sin^'d^  +  sin^cos^’d^ 
Therefore 

dl  _  cos^sin^  .  dip  sin^cos^ .  dip 
dh  sinl  dh  sinl  dh 


dl  _  cos#s imp  .  dip  sin^costftsinV' 
dh  -  sinl  dh  sinltan7 ( 1+h) 


Substituting  for  d^/dh,  Eq  (3.39),  gives 


-h/e 


dl  _  _L ,  cosasinqBe  '  _  cosqcosnfcost/’tan^ 

dh  Cp  ecos7sin7  (l+h)sin7 


+  2  u 


r  (l+h)ll/2 


ue. 


^^^(cos^tan7sin^  -  sin^) 


(3.4 


2 

2  (1+h)  .  ^  ^  ,  sintficoslsintf  ,0  ,, 

-  u  r  • — k — i — z — cosasinmcosmcostp  +  - — utt -  • — t  (3.4 

*  ug  tan7  r  r  r  cos0 (1+h) tan7sinl  v 


igure  11.  The  Osculating  and  Equatorial  Planes 
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The  last  term  cancels  with  the  second  term  and  with 


application  of  some  of  the  spherical  trigonometric 


equations,  Eqs  (3.42)  -  (3.47),  the  dl/dh  equation  is 


transformed  from  a  function  of  6,  <p ,  and  ip ,  to  a  function  of 


a  ,  fi ,  and  I . 


dl  L  cosasinaBe 


-h/e  fr  (l+h)l,/2 


dh  (L.  ecos7sin  7 


+  2  ui 


cosasinl , 
tan7 


(3.50) 


•  (cosat anor  -  sina)  -  w2  r  * — ^  — (sinlcoslcosasina) 
v  1  '  *  ug^tan7v  ' 


To  find  an  expression  for  da/dh  in  terms  of  the  desired 


variables,  Eq  (3.42)  is  differentiated. 


d <f>  _  coslsina  .  dl  sinlcosq  _  da 
dh  ~  cos <p  dh  cos^  dh 


Substituting  for  d<p/dh  with  Eq  (3.33)  produces 


da  _  _ 1 _  _  tana .  dl 

dh  ~  (l+h)tan7  tanl  dh 


(3.51) 


da  _  _ 1 _  _  ^L.  sinasinaBe 

dh  (l+h)tan7  ecos7sin7tanI 


-  2  w 


r^(l+h) 


sinacosl 


tan7 


(cosatan7  -  sina) 


2  (1+h)  .  2  2  -r 

+  u  r  — J — sin  a  cos  I 
*  ug  tan7 


(3.52) 


To  find  an  expression  for  dfl/dh  in  terms  of  the  desired 


variables,  Eq  (3.46)  is  differentiated. 


sin^  =  sinlcos(6-0) 


By  substituting  in  Eqs  (3.32)  and  (3.45),  this 


expression  can  be  rewritten  as 


d#  _  cosltan # .  dfi  _  cosltan#  coslcosq  .  dl 

dh  cos  ip  dh  (1+h) tan^cos#  cos^cos^  dh 


(3.53) 


Rewriting  Eq  (3.48)  gives  another  expression  for  dip/ dh : 


d#  _  sinl  .dl  _  sin#cos# 
dh  -  cos^sin ip  dh  (1+n)  tan^ 


(3.54) 


Equating  Eq  (3.53)  and  Eq  (3.54)  gives  dfl/dh  as  a  function 
of  dl/dh . 

dfi  _  cos  ip  cosltan#  _  sin#cos#cos# 

dh  cosltan#  (1+h) tan^cos#  (l+h)tan,y 


dl 

r  sinl 

coslcosal 

+  dh 

Lcos#sin# 

cos#cos#J 

dfi  _  _ cos# _ _  cosltanc 

dh  cos!tan#tan# (1+h) tan#  cos# 


coslsin^ 

2  * 
cos  # 


dl .  cos  ip  r< 
dh  cosIsin#L< 


s#cosaJ 


11  T  •  1  VWOU/N-WDW 

dh  coslsin^ Lcosa  r 

This  simplifies  to  the  following 

dfi  _  tana . dl 
dh  sinl  dh 


(3.55) 


Eq  (3.55)  can  be  "de-simplif ied"  by  substituting  Eq  (3.50) 
into  this  equation. 


,n  CT  n  -h/e  .  q  T r  (l+h)ll/2 

dfi  _  L.Be  '  sinasma  2w  *  v  ' 

dh  ~  Cp  esin,)(cos7sinI  tan7  ug^ 


•  [^sinacosatan^  -  sin  aj  -  ug  tan 7 — 'coslsin  a 


(3.56) 


The  du/dh  equation,  Eq  (3.34),  can  easily  be  transformed 
to  the  desired  orbital  variables  by  application  of  the 
spherical  trigonometric  relations. 


u 

2uB 

-We. 

(l+h) 

esin7 

ur* (1+h) 

1  /  2 

4wcosI 

— 

g* 

jsatan'yj 


2  2 

2  (l+h)  cosasin  Isina 

j  r  '  -» - L - - - 

*  g  tan7 


(3.57) 


The  last  equation  of  motion  to  be  transformed  is  the 
dq/dh  equation,  Eq  (3.38).  The  w2  term  of  Eq  (3.38)  is 

-  J2  r^q3  ^cos^sin^sin^tan7  +  cos2  <p j 

This  term  can  be  rewritten  as 

2  3  (l+h)2  .  .  2  Tj  .  2  T  •  2  1 

-  w  r  q  cosasinasin  Itan7  +  1  -  sm  Ism  a 

*  ug*  '  J 

The  dq/dh  equation  transformed  into  the  orbital  elements  is 
dq  q  fq2  ,1  CL  Be_h^E  „  [r» (1+h) 1 '  ^ 

dh  =  TuhTlu  "  Jj  “  — i — cosa  ~  2q“cosI[  Ug.  j  * 

-  w2  r^q3  +  ^  sinacosatan7sin2 1  +  1  -  sin2Isin2aj  (3.58) 


In  summary,  the  equations  of  motion  for  three- 
dimensional  rotating  planetary  entry  have  been  derived  for 
the  independent  variable,  non-dimensional  altitude,  and  some 


V 


convenient  dependent  variables  including  orbital 
inclination,  longitude  of  the  ascending  node,  and  argument 
of  latitude  at  epoch.  These  equations  are  given  below. 

du  -  _  _ H _ 2Bu  .  -h/ef,  Si.  cosataml 

dh  "  (1+h)  esin7  e  [  CD  cosatan7j 


-  4wcosI 


ur^ (1+h) 


2r  .  (1+h)2  cososin2  Isinq  f3 
*  g  tan7 


da _ a _ [a2_  il  _ 

dh  "  (1+h) Lu  XJ  cD 


CL  Be  h/e  „  T[r*(1 

Cn  e  ug 


r  (l+h)l,/2 


2  3  (1  + 

-  w  r  a  -*■ - 

*  ug 


dl  _  L . cosasingBe 
dh  -  CL.  ecos7sin7 


-h/e 


+  2ui 


r  (1+h)  11/2 


cosasinl 

tan7 


(cosatan7  -  sina) 


.2  _  .  (1+h) 


-  u  r  • — ^ — r — — (sinlcoslcosasina) 
*  ug  tan7  v 


CT  n  -h/e  .  .  „  f  r  (1  +  h)  1  ^2 

L.Be  '  si nasi no  2u;  * 

Cjj  esin7cos7sinI  tan7  ug^ 


sinacosatan7  -  sin2aj  - 


w2  r^ (1+h) 2 

— — t - *  cos  Is  in  a 

ug  tan7 


i  C,  .  .  n  -h/e 

_ 1 _  _  L.sinasingBe 

(l+h)tan7  Cp.  ecos7sin7tanI 


rJ1+h)  1 1  /2  sinacosl ,  .  .  . 

i  -  — r - (cosatan7-sina) 


XVvV.vW? 


[sinacosatan7sin2 I  +  1  -  sin2Isin2aj  (3. 


(3.6 


(3.6 


in  a cos  I 


(3.63) 


+  u  r 


*  ug  tan7 


Singularities  exist  for  these  equations  of  motion  for 
flight  path  angle  values  of  0.0  and  90.0  degrees  and  for  an 
orbital  inclination  angle  of  0.0  degrees.  Special 
consideration  must  be  taken  when  evaluating  Eqs  (3.59)  - 


The  rotating  Earth  terms  in  each  of  the  equations  of 
motion  are  examined  in  this  section.  This  examination  is 
conducted  for  orbital  inclination  angles  ranging  from  0.5  to 
75.0  degrees  and  vehicle  speeds  ranging  from  circular 
orbital  velocity  to  low  supersonic  speeds,  where  terminal 
maneuvers  such  as  landing  approaches  are  usually  initiated. 
The  rotating  Earth  terms  are  set  equal  to  zero  for  each 
equation  of  motion  and  are  then  checked  for  the  existance  of 
real  solutions.  The  existance  of  real  solutions  for  any  of 
these  expressions  indicates  trajectory  states  (specific 
combinations  of  values  of  u,  h,  q,  I,  0,  and  a  along  an 
entry  trajectory)  exist  where  that  particular  non-rotating 
Earth  equation  of  motion  is  valid.  In  a  later  section, 
these  trajectory  states  will  be  examined  in  more  detail. 

The  non-existance  of  real  solutions  to  the  rotating  Earth 
terms  in  any  one  of  the  equations  of  motion  indicates  the 
corresponding  non-rotating  equation  of  motion  will  be 
invalid  for  any  Earth  entry  trajectory.  A  solution  that 
accounts  for  the  Earth’s  rotation  is  needed  for  any  non¬ 
rotating  Earth  equation  of  motion  that  is  invalid  for  the 
full  range  of  inclination  angle  and  speeds  investigated. 
These  solutions  are  developed  in  Section  V  by  application  of 
the  method  of  matched  asymptotic  expansions. 


The  Second  Small  Parameter 


To  begin,  the  equations  of  motion  are  placed  in  a  more 
convenient  form  to  work  with.  In  Section  III,  £  was  shown 
to  be  approximately  equal  to  the  ratio  of  the  atmospheric 
scale  height  and  the  radius  of  the  planet.  A  convenient 
second  small  parameter  is  introduced  to  replace  the  w  terms 
in  the  equations  of  motion.  This  second  small  parameter  is 
defined  as  the  square  of  the  ratio  of  the  planet’s 
rotational  velocity  at  the  equator  and  the  prograde, 
equatorial  circular  orbital  velocity  at  the  surface.  The 
orginal  small  parameter,  e,  is  re-labelled  e ^  and  the  second 

small  parameter  is  labelled  e . 


-  €2  (1  +  h)  (tan^sinacosasin  I  +  1  -  sin  Isin  a)  (4.3) 


CT  D  -h/e 

L  Be  l  sinasma 


da  1 _ L.  Be 

dh  (1+h) tan 7  e  ^ 


cos'ysin'ytanl 


1/2 

n\  Q+h)l  sinacosl/  ,  .  % 

-  2  e  A - *-  - t - (cosatani  -  sina) 

.2  u  tan7  v  '  ' 


+  (  h);.sin!acos!l 

2  utan7 


(4.4) 


=  gL.Be~h/ei  singsina  +  2f  ii±hll 1  X! 5isa (cOBat  -sina) 

dh  Cp  e^  cos'jsin'jsmi  L  2  u  tan7 x  ' 


2 

x 4  ,v2 . sin  acosl 
—  e„ (l+h)  7 

2  utan^ 


,T  CT  D  -h/e 
dl  L  Be  '  l  sxnacosa 


(4.5) 


dh  cD  Cj 


cos^sin^ 


L  2  u 


l/2 

+h)  1  '  . cosasinl , 

u  J  tan^ 


x  ,  .  \  „  m2  cosasinacoslsinl 

•  (cosatan'y  -  sina)  -  C2  (1+h)  - - -  - 


(4.6) 


It  was  previously  noted  that  the  latter  three  of  the 


above  equations  are  coupled 


dfi  _  tang  .  dj[ 
dh  sinl  dh 


da  _  _ 1 _  _  tana  .  dl 

dh  ~~  (l+h)tan7  tanl  dh 


The  Rotating  Earth  Terms 


(4.7) 


(4.8) 


The  equations  of  motion  for  lifting  atmospheric  entry 
for  a  rotating  Earth  can  be  written  as  the  sum  of  the 
equations  for  a  non-rotating  Earth  and  of  the  terms  that 


Nonrotate 


(1+h)  Lu  AJ  C 


D  ei 


•  cosct 


=  -  COSl 


1  Rotat 


(4.10) 


2  Q  2  2  2 

-  e  (1  +  h)  (tan-ysinacosasin  I  +  1  -  sin  Isin  a) 


Argument  of  Latitude  at  Epoch.  The  rotating  Earth 
equation  of  motion  for  argument  of  latitude  at  epoch,  given 
by  Eq  (4.4),  can  be  rewritten  as 


=  Nonrotate  +  Rotate 


where 


da  =  1 _ L.Be 

dh  Nonrotate  (l  +  h)t.an7  “  Cp  £ 


CT  n  -h/e 

L  Be  '  l  sinasina 


cos7sin7tanI 


I  Rotat' 


1/2 

« r  (l+h)l  '  sinacosl,  ,  .  x 

=  -  2[£2-i-^]  tan7  “  ■»■“) 


*  c  (l^h)?-sin?°cos21 
2  utan7 


(4.11) 


Longitude  of  the  Ascending  Node.  The  rotating  Earth 
equation  of  motion  for  longitude  of  the  ascending  node, 
given  by  Eq  (4.5),  can  be  rewritten  as 


=  Nonrotate  +  Rotate 


where 


dn 

dh 


Nonrotate 


L.Be  '  l  sing sing 
Cp  cos7sin,ysinI 


dfi 

dh 


Rotate 


.sing 
tan  7 


(cosgtan7 


sing) 


,,  .IX2  . sin2  gcosl 
£  ( 1 +n  )  t 

2  utan7 


(4.12) 


Inclination  Angle  Equation.  The  rotating  Earth  equation 
of  motion  for  orbital  inclination  angle,  given  by  Eq  (4.6), 
can  be  rewritten  as 

=  Nonrotate  +  Rotate 


where 


dl 

dh 


CL  Be~h/£l  .si _ 

Nonrotate  CD  £t  cos7sin7 


nacosg 


dl 

dh 


-  2 

[e  -^11 

Rotate 

.2  u 

tan7 


(cosgtan7  -  sing) 


,  n2  cosgsingcoslsml 

e  (1  +  h) - t - 

2  utan7 


(4.13) 


Investigation  of  Real  Solutions  to  the  Rotate  Terms 

Eqs  (4.9)  -  (4.13)  are  the  Rotate  equations,  those  parts 


of  the  complete  equations  of  motion  which  account  for  the 
rotating  Earth.  In  the  following  pages,  the  five  Rotate 
equations  are  each  set  equal  to  zero  and  then  examined  for 


1 


1 


From  Eq  (4.9) 


Rotate 


-  4^e2u(l+h)j  cosl 


2e2  (1+h) 


2 

2  cosgsin  Isina 
tan 7 


Setting  this  equation  equal  to  zero  and  solving  for  one  of 
the  five  variables  gives 


7  =  tan 

Rotate  =  0 


2 

-1  -Ccososinqsin  I 
2cosI 


(4.1 


where  c=p^%J,/2 


(4.1 


Substitution  of  a  few  realistic  combinations  of  u,  h, 
and  a  indicate  real  solutions  exist  for  Eq  (4.14). 

Rotate  Term  Solutions  for  the  Inclination  Ancle 


Equation .  From  Eq  (4.13) 


3K  _ _ =  z[e2X1?1]  '“t”'”1!'05111"’  “  sin“> 


Rotate 


-  e2  (1+h) 


2  . cosa si  no cos Is ini 
utan^ 


Setting  this  equation  equal  to  zero  and  substituting  in  Eq 
(4.15)  for  C  gives 


O  *  T  S  X  nO  V  .  .  -y  -r 

zsinicosa  cosa  -  .  -  =  t - sinacosasinlcosl 

tan7J  tan7 


(4.16) 


One  set  of  trivial  solutions  to  this  equation  is  in  the  form 


sinlcosa  =  0  for  any  C  and  7 


k». 


The  w  term  contribution  to  the  dl/dh  equation  of  motion  is 
therefore  equal  to  zero  when  sinl  and/or  cosa  =  0.  However 
due  to  the  fact  the  equations  of  motion  are  singular  for 
sinl  =  0  as  discussed  in  Section  III,  the  trivial  solution 
is  actually 

cos  a  =  0  (a  =  .  .  .  -3ir/2,  -n/2,  r/2,  3tt/2,  ...) 

Nontrivial  solutions  to  Eq  (4.16)  can  be  found  by  solving 
for  one  of  the  five  variables. 


7  =  tan 

Rotate  =  0 


l  sing  +  Ccoslsina 
2cosg 


(4.17 


Substitution  of  a  few  realistic  combinations  of  u,  h,  I,  an 
a  indicate  real  solutions  exist  for  Eq  (4.17). 

Rotate  Term  Solutions  for  the  Argument  of  Latitude  at 


Epoch .  From  Eq  (4.11) 


I  Rotate 


l/2 

2  (1+h) '  '  . sing  cosl 

^12  u  tan^ 


(cosgtan'j  -  sing) 


2  9 

, ,  L\2  sin  g  cos  I 

+  e  (1+h) - 7 - 

2  utan^y 

Setting  this  equation  equal  to  zero  and  substituting  in  Eq 
(4.15)  for  C  gives 


o  •  t  sing] 

2smgcosl  cosg  -  ^an~ j  = 


tan7J  tan^ 


C  .2  2  T 

- -‘sin  geos  1 


(4.18) 


One  set  of  trivial  solutions  to  this  equation  is  in  the  form 


sinacosl  =  0  for  any  C  and  7 


The  u  term  contribution  to  the  da/dh  equation  of  motion  is 
hence  equal  to  zero  when  sina  and/or  cosl  are  zero. 

However,  due  to  the  range  of  inclination  angle  selected  for 
study,  cosl  is  never  equal  to  zero  and  the  trivial  solution 


sina  =  0  (a  =  .  .  .  -2?r,  -ir,  0,  ir ,  2ir ,  •••) 

Nontrivial  solutions  to  Eq  (4.18)  can  be  found  by  solving 
for  one  of  the  five  variables  in  this  equation. 


7  =  tan 

Rotate  =  0 


-1  Ccoslsinq  +  sing 
2cosa 


(4.19) 


Substitution  of  a  few  realistic  combinations  of  u,  h,  I,  and 
a  indicate  real  solutions  exist  for  Eq  (4.19). 

Rotate  Term  Solutions  for  the  Longitude  of  the 
Ascending  Node.  From  Eq  (4.12) 


=  2 


1  /  2 


Rotate 


(1+h)!'  sina,  ^  N 

[e2^rj  t^,(cosat*n'i  ~  sinQ> 


e2  (l+h) 


■  2  T 
2  .sin  a  cosl 

utan7 


Setting  this  equation  equal  to  zero  and  substituting  in  Eq 
(4.15)  for  C  gives 


0  .  sina  n  ■  2  cosl 

2sina  cosa  -  t -  =  Csin  a-— — - 

tan7_  tan7 


(4.20) 


One  set  of  trivial  solutions  to  this  equation  is  in  the  form 
sina  =  0  for  any  C,  I,  and  7 


Hence,  the  u  term  contribution  to  the  da/dh  equation  of 
motion  is  equal  to  zero  for  a  =  ...  -2ir ,  -* ,  0,  if,  2 k, 


Nontrivial  solutions  to  Eq  (4.20)  are  in  the  form 


-l 

Ccoslsina  +  sina 

7 

=  tan 

Rotate  =  0 

2cosa 

(4.21) 


Substitution  of  a  few  realistic  combinations  of  u,  h,  I,  and 
a  indicate  real  solutions  exist  for  Eq  (4.21). 

Coupling  of  Three  Rotate  Term  Solutions.  The  non¬ 
trivial  solutions  for  the  rotating  term  expressions  are 
identical  for  the  da/dh,  df]/dh,  and  dl/dh  equations,  Eqs 
(4.17),  (4.19),  and  (4.21).  This  can  also  be  seen  by  noting 

that  the  dl/dh,  dfl/dh,  and  da/dh  equations  are  coupled  by 
Eqs  (4.7)  and  (4.8). 


dfi  _  tana . dl 
dh  sinl  dh 


(4.7) 


da  _  _ 1 _  _  tana  .  dl 

dh  (l+h)tan7  tanl  dh 


(4.8) 


These  coupling  relations  can  be  written  in  terms  of  the 
Nonrotate  and  Rotate  expressions: 


dfi 

dh  Nonrotate 
+  Rotate 


tana . dl 
sinl  dh 


+ 

Nonrotate 


tana . dl 
sinl  dh 


Rotate 


(4.22) 


da 

dh  Nonrotate 
+  Rotate 


1 

(1 +h) tan^ 


tana , dl 
tanl  dh 


Nonrotate 


tana . dl 
tanl  dh 


Rotate 


(4.23) 


Since  the  Rotate  terms  contain  u  or  e by  definition,  the 

first  term  of  Eq  (4.23)  must  be  a  Nonrotate  expression. 
Hence , 


da 

dh 

In 


Rotate 

addition 


tang . dl 
tanl  dh 


Rotate 


(4.24) 


dn 

dh 


Rotate 


tang . dl 
sinl  dh 


Rotate 


(4.25) 


From  Eqs  (4.7)  and  (4.8)  it  can  be  seen  that  the  dfl/dh 
and  dg/dh  equations  differ  significantly.  However, 
comparing  like  terms  indicates 


da 

dh 


Rotate 


cos 


Rotate 


(4.26) 


These  results  simplify  the  search  for  the  trajectory 
states  where  the  non-rotating  Earth  equations  of  motion  are 
valid.  Instead  of  detailed  examination  of  the  non-trivial 
solutions  of  all  three  of  the  da/dh,  dQ/dh,  and  dl/dh  Rotate 
equations,  examination  of  only  one  of  them  is  required.  The 
trivial  solutions  for  rotating  term  expressions  for  the 
dfl/dh  and  da/dh  Rotate  equations  are  the  same,  sina  =  0,  and 
the  trivial  solution  for  the  rotating  term  expression  for 


the  dl/dh  Rotate  equation  is  cosa  =  0. 

Rotate  Term  Solutions  for  the  Flight  Path  Angle 


dh 


Rotate 


COS1 


o 

2  Q  •  2  2  2 

-  e2  (1+h)  (tan7sinacosasin  I  +  1  -  sin  Isin  a) 


Setting  this  equation  equal  to  zero  gives 


3  2  2  2  1 

0  =  2qA^ cosl  +q  tan7sinacosasin  I  +  1  -  sin  Isin  al  (4.27 


where 


A  = 
1 


T  1  /2 


e2  (1+h)' 


(4.28 


Because  singularities  exist  for  I  or  7  =  0.0  or  90.0 

degrees,  q  =  0  can  not  be  considered  as  a  trivial 

solution.  To  examine  possible  solutions  to  Eq  (4.27),  some 

additional  variables,  A.,  are  defined  for  convenience. 

’  x 


let 

A2 

1  .  2  T  .  2 

=  1  -  sin  Isin  a 

(4.29 

A3 

.  2 

=  sinacosasin  I 

(4.30 

A 

A 

=  2A^  cosl 

(4.31 

Using  these  variables,  dividing  by  q,  and  remembering 
q  =  COS7,  equation  Eq  (4.27)  can  be  rewritten 

A  +  q2  A  ta.n'j  +  q2  A  =0 

*1  u  Z 

2 

A  +  A  cos  7  =  -  A  sin7cos7  (4.32 

H  Z  o 


Squaring  this  expression 


2 


A  2  cos4  7  +  A  +  2A  A  cos2  7  -  A  2  cos2  7  +  A  2  cos4  7  =  0 
2  *  4  4  2  '  3  '3  ' 


(A  2  +  A  2)i/2  +  (2A  A  -  A  2)i/  +  A  2  =  0 

2  3  42  3  4 


(4.33 


,  2  2 

where  1/  =  q  -  cos  7 


(4.34 


Eq  (4.33)  can  be  easily  solved  in  the  form 


Ax/  +  Biv  +  C  =  0 


where  A  =  A?  2  +  A3  2  ,  B  =  2A^  A?  -  A^  ,  C  =  A^  2  ,  and 


-B  *  [b2  -  4Acl 
"  =  - 2A - 

The  Rotate  term  solution  for  the  flight  path  angle  equation 
is  hence  given  by 


-2A  A  +  A  2  *  [a. 


-  4A  A  A/  -  4A.  A 


At  this  point  1/  has  been  placed  in  terms  of  the 
variables  A^  ,  A^  ,  A^  ,  and  A^  .  Looking  at  each  of  these 

expressions  individually  aids  in  the  determination  of 
general  trends  of  the  values  of  u . 

A  Equation.  It  can  be  easily  seen  that  A  has  a 

maximum  value  of  1.0  when  sina  =  0,  for  any  value  of  sinl 

within  the  inclination  range  of  interest,  0.5  I  75.0 

degrees.  A  minimum  of  0.0670  occurs  for  A  when  I  =  75.0 


degrees  and  a  =  90.0  degrees  or  90.0  ±  180.0  degrees. 
Thus,  A2  is  never  negative. 


A  Equation.  A  has  a  maximum  value  of  0.466  for 

—3 -  -  3 

I  =  75.0  degrees  and  a  =  45.0  degrees.  A  minimum  value  of 


-0.466  occurs  when  I  =  75.0  degrees  and  a  =  -45.0  degrees 


A^  Equation .  Since  A^  is  always  positive,  it  can 


be  easily  seen  A  has  a  maximum  value  of  34.0  for 

4 


I  =  0.5  degrees,  u  =  1.0,  and  h  =  0.0  .  A  has  a  minimum  of 

4 


0.0  for  u  =  0.0  .  Minimum  values  of  A  for  non-zero  u  occur 

4 


for  I  =  75  degrees.  For  a  prograde  orbit  A  is  always 

4 


positive 


Solutions .  For  real  solutions  to  exist  in  Eq 
(4.35),  and  hence  for  the  rotating  terms  in  the  dq/dh 
equation  to  ever  have  a  zero  contribution,  the  expression 
within  the  square  root  must  be  positive.  Factoring  A3  out 

of  the  root  expression  in  Eq  (4.35)  leaves  the  following 
condition  for  the  existence  of  real  roots 


A  2  -  4 [ A  A  +  A  21  >  0 

3  L  2  4  4  J 

where  A?  and  A^  are  always  positive  and  A^ 


(4.36) 


can  be  positive 


or  negative  over  the  range  of  inclination  angle  examined. 


Vvvv\ v\ 


r«* 

fw 

JV 

K 


Eq  (4.36)  indicates  real  solutions  to  the  dq/dh  Rotate 
equation  do  not  exist  for  realistic  entry  trajectories 
except  for  values  of  u  corresponding  to  low  speed,  subsonic 
and  near  subsonic  flight.  Therefore,  no  real  roots  to  Eq 
(4.35)  exist  for  orbital  inclination  angles  between  0.5  and 
75.0  degrees  and  vehicle  speeds  ranging  from  circular 
orbital  velocities  to  low  supersonic  speeds.  This  result 
indicates  the  non-rotating  Earth  dq/dh  equation  of  motion  is 
invalid  for  Earth  atmospheric  entry  for  the  ranges  of 
orbital  inclination  angle  and  velocity  investigated. 
Therefore,  a  first  order  solution  (Section  V)  to  the 
complete  dq/dh  equation  is  required  for  Earth  atmospheric 
entry  analysis.  This  solution  is  developed  in  the  next 
section  using  the  method  of  matched  asymptotic  expansions. 
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Solutions  to  the  Equations  of  Motion 


Matched  Asymptotic  ExDansions 


In  Section  IV,  the  terms  in  the  five  equations  of  motion 


that  account  for  Earth’s  rotation  were  examined.  It  was 


determined  that  trajectory  states  exist  for  four  of  the  five 


equations  of  motion  where  the  rotating  Earth  terms  give  a 


zero  contribution.  For  the  fifth  equation  of  motion,  the 


dq/dh  equation,  it  was  shown  that  the  rotating  earth  terms 


always  have  a  non-zero  contribution  for  realistic  lifting 


entry  trajectories  with  entry  inclination  angles  between  0.5 


and  75.0  degrees  and  speeds  ranging  from  circular  orbital 


velocities  to  low  supersonic  speeds.  Therefore,  to 


adequately  model  Earth  atmospheric  entry,  a  solution  to  the 


dq/dh  equation  is  required  that  includes  the  rotating  Earth 


effects.  This  solution,  along  with  the  non-rotating  Earth 


solutions  to  the  other  four  equations  of  motion,  is 


developed  in  this  section,  by  application  of  the  method  of 


matched  asymptotic  expansions. 


Combining  the  Small  Parameters 


In  Section  III,  e  was  defined  to  be  the  ratio  of  the 


planetary  atmosphere  scale  height  and  the  planet’s  radius. 


A  second  small  parameter  was  introduced  in  Section  IV  to 


non-dimensional i ze  the  rotating  planet  terms  in  the 


equations  of  motion.  This  second  small  parameter  was 


defined  as  the  square  of  the  quotient  of  the  planet’s 


rotational  velocity  at  the  equator  and  the  prograde, 


equatorial  circular  orbital  velocity  at  the  mean  planetary 
radius . 


e .  =  a-nd  e„  = 

l  p  r  2 


The  general  equations  of  motion  for  three-dimensional, 
rotating  planetary  entry  were  given  in  Section  IV  in  the 


(l+h) 


2uBe  h/£i  .  F. 
sin7 


1  tan^cosaj 


l/2  2 

x  n.  1  \2cosasin  Isina 


-  4  e2u(l+h)j  cosl  -  2e2(l+h)' 


tan7 


(5.2) 


‘  Ti^hytu”  -  x]  -  ' COBl 

-  e  (l+h)2 •^-(tan7sinacosasin2 I  +  1  -  sin2Isin2a)  (5.3) 


da  =  1 _ L.  Be 

dh  (l+h)tan7  Cp 


Ct  n  -h/c, 

L  Be  1  sinasma 


cos7sin7tanI 


:  f£  Xi±hll 

2  U 


sinqcosl 

tan7 


(cosatan7  -  sina) 


+  €  (l+h)8-5in?°C°S;i 

2  utan7 


(5.4) 


dfl  _  ^L.Be  sina 

dh  Cp  e)  cos7sin7 


in^sinl  +  2[£2-Q^lli]  /'|^(=°Mta„7-sina) 


dl  _  .  Be  singcosa  +  ^  [ e  (l+h)1  1  cosasina . 

dh  _  Cp  e ]  cos7sin7  ^L£2  u  J  tan7 

•  (cosatan7-sina)  -  (1+h) - utan7 -  (5-6) 

Values  for  the  two  small  parameters  must  be  calculated 
and  substituted  into  Eqs  (5.1)  -  (5.6)  for  each  planet 
studied.  For  the  Earth,  the  value  of  the  first  small 
parameter  is  approximately  1/900  and  the  value  of  the  second 
small  parameter  is  approximately  1/289  .  To  simplify  the 
following  analysis,  it  is  noted  that  the  first  small 
parameter  is  approximately  equal  to  one  third  of  the  second 
small  parameter.  Hence,  a  new  small  parameter,  e,  is 
defined  for  Earth  entry  analysis  as 

e  =  =  e,,/2=  <v3)1/2;  ei  =  g2  and  s  =  3e2  (5’7) 


The  Equations  of  Motion  for  Earth  Atmospheric  Entry .  To 
obtain  the  equations  governing  Earth  atmospheric  entry,  the 
first  and  second  small  parameters  are  replaced  in  the 
equations  of  motion  by  e  and  a  constant.  The  resulting 
equations  are  used  in  the  derivation  of  solutions  for  Earth 
entry : 


iu  -  u 


2uBe 


dh  “  (1+h)  2  . 

v  '  e  sin7 


-h/ e  r  Cj  1  r  i 1  /2 

-•  ^1  +  ^=tan7cosaJ  -  4e^3u(l+h)J  cosl 


D 


.  2  T  • 

e  2  /,  ,,_\2  .cosasin  Isina 
-  6e  (1+h) - t - 


(5.8 


MJ 


dh  "  (l  +  h) Lu  AJ  CD  e2 


coso  -  - - -j 


O 

2  2  Q  •  2  2  2 

-  3e  (l+h)  •^-(tan'ysinacosQsin  I  +1-  sin  Isin  a)  (5.C 


,  .  CT  D  -h/e2 

da  _  1 _  _  L. Be _ singsina 

dh  (l  +  h)tan7  £2  cos-ysii^tanl 

0  f3(l+h)l 1 sina  cosl,  .  .  x 

_  2e  [“  ^ - L\  - — (cosatan7  -  sina) 

+  ,  (  h)g. sin2 a  cos2I 

v  '  utan^ 


(5. 1C 


-h/c' 


dn  L  Be  ^  sinasma  n  T3(l+h)  ^2sina,  ,  . 

jr  =  rr~  — ; - : - ; — =f  +  *£  I — ^ -  '  t - ( cosatani-s  i  na  ) 

dh  Cp  2  cos7s1n7s1.nl  L  u  tan^ v  ' 


+  3£2(i+h)2  si°’;“sl 


(5.  n 


,T  Cj  r,  -h/e  . 
dl  _  L . Be  '  _ singcosg 

dh  Cp  ^2  cos^sin 7 


1  /2 

ngcosg  +  T3(l+h) 1  '  . cosasina 

sTfsimr  e  L  u  J  tan7 


✓  ,  •  \  -a  2  ,  v2  cosasinacoslsinl 

•  (cosatan7  -  sma)  -  3e  (l+h) - r - 

u  u2l  n  j 


(5.12 


It  was  previously  seen  that  the  latter  three  of  the 


above  equations  are  coupled: 

dn  _  tang . dl 
dh  sinl  dh 


da  _  _ 1 _  _  tana  .  dl 

dh  (l+h)tan7  tanl  dh 


(5.12 


(5.14 


Eqs  (5.8)  -  (5.12)  describe  the  three-dimensional, 
rotating  planetary  entry  of  Earth.  However,  the  generi< 
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application  of  these  equations  to  some  of  the  other  planets 
is  not  lost.  As  an  example,  consider  the  problem  of  Mars 
atmospheric  entry.  To  make  Eqs  (5.1)  -  (5.6)  applicable  for 
Mars  entry  analysis,  the  first  and  second  small  parameters 
are  calculated  using  Mars  planetary  constants.  A  new  small 
parameter,  e,  differing  from  Earth’s  by  only  a  constant  is 
then  defined.  For  Mars,  the  first  small  parameter  has  a 
value  of  approximately  1/350  .  The  second  small  parameter 
is  calculated  from  IAU  defined  constants  to  be  approximately 
1/218. 

11  2r*  _  1 
£1  "  /Jr  ~  350  e2  "  W  g^  “  218 

The  first  small  parameter  is  approximately  equal  to 
three  fifths  of  the  second  small  parameter.  Therefore,  a 
unique  e  is  defined  for  Mars  entry  analysis: 

e  =  1^7  =  (e,)l/2=  (3e,  /5) 


or 


e2  and  e 

2 


(5.15) 


To  adjust  the  equations  of  motion  for  Martian 
atmospheric  entry,  the  new  first  and  second  small 
parameters,  as  a  function  of  e  and  a  constant,  are 
substituted  into  Eqs  (5.1)  -  (5.6).  Hence,  the  equations  of 
motion  for  Mars  entry  will  differ  from  those  for  Earth  entry 
analysis  by  only  a  constant  in  each  expression  containing  e . 
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The  Method  of  Matched  Asymptotic  Expansions 


While  many  analytical  methods  have  been  applied  to  the 
problem  of  planetary  atmospheric  entry,  one  method  remains 
relatively  unexploited.  This  is  the  method  of  matched 
asymptotic  expansions.  From  an  entry  vehicle’s  viewpoint, 
planet’s  atmosphere  forms  a  boundary  layer  of  density  in 
space.  Aerodynamic  forces  on  the  vehicle  change  from 
insignificant  to  dominating  for  a  relatively  small  change  i 
altitude.  A  variety  of  methods,  including  composite 
expansions,  multiple  scales,  and  matched  asymptotic 
expansions,  have  been  utilized  to  solve  boundary  value 
problems.  However,  the  method  of  matched  asymptotic 
expansions  is  more  "versitle"  and  "effective"  than  these 
other  methods  for  both  linear  and  nonlinear  problems 
composed  of  partial  or  ordinary  differential  equations 
(Nayfeh,  1985:257,258).  Past  papers  describe  successful 
applications  of  matched  asymptotic  expansions  to  solve 
limited  flight  mechanics  problems  involving  lifting 
atmospheric  entry.  (Busemann  and  others: 1976,  Shi  and 
Pottsepp : 1969 ,  Shi: 1971,  Shi  and  others: 1971,  Willes  and 
others : 1967) .  Most  of  these  efforts  used  two-dimensional 
equations  of  motion  and  assumed  a  spherical ,  non-rotating 
Earth  model  except  Busemann  (Busemann  and  others,  1976), 
where  the  three-dimensional  equations  of  motion  were 
employed.  In  the  following  pages,  the  solution  for  the 
dq/dh  equation  of  motion  for  a  rotating  Earth  is  derived. 

To  obtain  this  solution,  the  solutions  for  the  non-rotating 


equations  of  motion  are  first  derived  from  the  rotating 
Earth  equations  of  motion. 

For  a  boundary  layer  problem,  the  method  of  matched 
asymptotic  expansions  gives  two  or  more  solutions,  each 
valid  in  specific  regions  of  the  domain.  These  solutions 
have  some  overlap  and  can  therefore  be  matched.  The 
matching  conditions  allow  for  the  construction  of  a 
composite  solution  which  is  valid  over  the  entire  domain. 

In  this  analysis,  the  entry  problem  is  treated  as  a  boundary 
layer  problem  with  a  solution  developed  away  from  the 
boundary  layer,  valid  for  exo-atmospheric  flight,  and  a 
solution  developed  at  the  boundary  layer,  valid  for 
atmospheric  flight.  These  two  solutions  to  the  equations  of 
motion  can  be  independently  applied  in  their  respective 
domains.  However,  it  is  more  convenient  to  have  a  solution 
valid  for  the  full  range  of  entry  conditions.  Hence,  the 
inner  and  outer  solutions  are  matched  and  a  composite 
solution  valid  in  both  regimes  is  constructed. 


Outer  Expansions 

The  solution  developed  for  the  exo-atmospheric  portion 


of  the  domain  is  called  the  outer  solution.  This  solution 


is  developed  from  asymptotic  expansions  of  the  equations  of 
motion  using  the  small  parameter,  e.  The  outer  solution 
variables  are  denoted  by  the  superscript  "o".  Straight¬ 
forward  expansions  of  the  outer  solution  variables  are  made 
and  are  as  follows: 


V\ 


* £  'l  ♦ £  "l  4 


%  +  £'q.  +  *'%  + 


I  +  e1 1  +  e2I  + 


(5.16) 


n  +  e‘n.  +  e2n„  + 


a  +  £ 1  a  +  e2a  + 

0  1  2 


nr  +  e1  <>  +  e2  nf  + 

'0  '1  '2 


The  solutions  for  lifting  atmospheric  entry  of  a  non- 


rotating  planet  are  of  order  £  .  Adding  the  rotating  planet 


model  to  the  atmospheric  entry  problem  causes  additional 


terms  of  order  £  and  £2  to  appear  in  the  expansions  of  the 


equations  of  motion.  These  additional  terms  are 


relatively  complex  but  require  the  solutions  to  be  carried 


out  to  order  £  in  order  to  account  for  Coriolis 


acceleration  on  the  flight  vehicle.  To  also  model  the 


usually  insignificant  Centripetal  or  Transport 


acceleration,  solutions  would  have  to  be  carried  out  to 


order  £  .  The  solutions  developed  in  this  study  are  to 


order  e°  for  Eqs  (5.8),  (5.10),  (5.11),  and  (5.12)  and  to 


order  £ 1  for  Eq  (5.9),  the  dq/dh  equation  of  motion. 


Solutions  developed  to  order  £J  act  as  a  correction  to  the 


zero  order  solutions,  a  correction  which  accounts  for  the 


significant  effects  of  a  rotating  planetary  atmosphere. 

In  the  following  pages,  the  equations  of  motion  are 
expanded  for  small  e  using  the  straight-forward  expansions 
given  by  Eq  (5.16).  Many  of  the  common  expressions  that 
contain  e  in  the  equations  of  motion  are  expanded  in  more 
detail  in  Appendix  B. 

The  du/dh  Outer  Expansion.  The  outer  expansion  of  the 
du/dh  equation  is  given  by 
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The  dl/dh  Outer  Expansion.  The  outer  expansion  of  the 
dl/dh  equation  is  as  follows: 
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The  dfl/dh  Outer  Expansion.  It.  has  been  noted  that,  the 
dO/dh  equation  is  a  function  of  the  dl/dh  equation  as  well 
as  a  function  of  a  and  I.  This  dependence  exists  for  the 
outer  expansions  of  these  equations  as  well.  From  Eq  (5.13) 
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The  outer  expansion  of  the  dfj/dh  equation  is 
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The  da/dh  Outer  Expansion.  It  has  been  noted  that  the 
da/dh  equation  is  a  function  of  the  dl/dh  equation  as  well 
as  a  function  of  a,  h,  7,  and  I.  This  dependence  exists  for 
the  outer  expansions  of  these  equations  as  well.  From 


Eq  (5.14) 
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The  outer  expansion  of  the  da/dh  equation  is 
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The  dq/dh  Outer  Expansion.  The  outer  expansion  of  the 
dq/dh  equation  is  as  follows: 
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Solutions  to  this  set  of  differential  equations  are 
derived  in  Appendix  C  and  are  given  below.  Eq  (5.30)  is  the 
outer  zero  order  e  solution  to  the  dq/dh  equation  of  motion. 
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e1  Terms.  The  complete  set  of  e1  term  outer  expansion 
differential  equations,  derived  above,  are 
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A  complete  solution  to  the  dq/dh  equation  requires  that 
the  solutions  to  both  Eqs  (5.18)  and  (5.26)  be  found. 
Solutions  to  these  differential  equations  are  derived  in 
Appendix  C  and  are  given  below.  Eqs  (5.32)  and  (5.33)  are 
the  solutions  to  the  order  e  terms  in  the  outer  solutions  of 
the  du/dh  and  dq/dh  equations,  respectively. 
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where  C  ,  C  ,  and  C  are  functions  of  C  and  C  ,  and  C 
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is  a  constant  of  integration.  Expressions  for  these 

constants  are  given  in  Appendix  C. 

Order  e  Outer  Solution  for  the  dq/dh  Equation  of  Motion. 
The  order  e  outer  solution  for  the  cosine  of  the  flight  path 
angle  is  given  by  the  substitution  of  Eqs  (5.30)  and  (5.33) 
into  Eq  (5 . 16) . 
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Inner  Expansions 

The  solution  developed  for  the  atmospheric  portion  of 
the  domain  is  called  the  inner  solution.  This  solution  is 
developed  from  asymptotic  expansions  of  the  equations  of 
motion  using  the  small  parameter,  e.  The  outer  solution 
variables  are  denoted  by  the  superscript  "i".  Straight¬ 
forward  expansions  of  the  inner  solution  variables  are  made 
and  are  given  below. 
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To  derive  the  inner  expansions,  a  new  independent 
variable  is  required  so  as  to  force  the  equations  of  motion 
to  focus  on  the  boundary  layer  region.  £,  a  magnified 
version  of  h,  is  the  inner  expansion  independent  variable. 
Hence,  £  is  defined  as  the  magnified  non-dimensional 
altitude  and  it  is  given  by 
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Therefore  ^  =  e2  and  by  the  chain  rule, 
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The  equations  of  motion  for  Earth  atmospheric  entry  are 
now  transformed  from  functions  of  h,  represented  by  Eqs 
(5.8)  -  (5.12),  to  functions  of  the  magnified  non- 
dimensional  altitude,  £. 
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The  latter  three  of  the  above  equations  are  coupled 
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In  the  following  pages,  these  equations  of  motion  are 
expanded  for  small  e  using  the  straight-forward  expansions 
of  Eq  (5.35) .  Many  of  the  common  expressions  that  contain  e 
are  expanded  in  detail  in  Appendix  B. 

The  du/d£  Equation  Expansion.  The  inner  expansion  of 
the  du/dh  equation  is  given  by 
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The  dl/d£  Equation  Expansion.  The  inner  expansion  of 
the  dl/d£  equation  is 

dl  dl  CL 

-3-7^-  +  ejr-  +  0(e2)  =  7^— ’Be  ^sina(cosa  -  eo  6ina  )• 

d£  Op.  0  1  0 


1  sin70  1  cos7c 

-  +  £7  -  — ; -  -  £7  - 

COS7  '1  2  sin7  '1.2 

'n  ms  *n  kid  <-* 


cos  7 


sm  7n 


,1/2 

2e3  3(1  +  e2£)  — - e-r]  (COSQ0  -  ea  sina  ) 

0  u 


(sina„  +  ea , cosa  )  — — 
v  0  1  0  J  L tan 


t - -  J  (cos“o  -  e“,ein °„> 

'0  sin  7„  L 


[tan7  +  - - - - 1  -  (sina  +  ea  cosa  ) 

sin  7 


-  3e 


‘“♦'’o’lK:  -  ^[ts 


n  2  J  ltan7„  •  2 

0  u  '0  sm  7„ 

0  '0 


(eina  +  ea,  cosaj  (cosa„  -  £a,  sinaj 


,1 


'  (sinIQ  +  elj  cosIQ )  (coslo  -  e^sin^)  +  0(e  ) 


This  equation  reduces  to 


dl 

0 

dl 

+ 

1 

COS 

sin70  i 

+  e'y  -  — = — ■ 

'l  2  _  smr) 
cos  7  .  'i 

In 


n 7^  '1.2  v  ' 

'0  sin  7  . 


dl  CT 


o  .  “■‘■o  L  n  -f  . 

e  terms:  j-t—  =  j=r-’Be  ^sintr1 

dC  CD 


6  terms : 


Be  ^*sina  -a 


l  - : -  +  7 

l  cos70  sin^  '  l 


n1o 

(5.47) 

(5.48) 

cos  a 

0 

COSO 

0 

cos2  7^ 
'0 

^1  .  2 
sin  7o 

- 

The  d0/d£  Equation  Expansion.  It  has  been  noted  that 
the  dfl/d£  equation  is  a  function  of  the  dl/d£  equation  as 
well  as  a  function  of  a  and  I.  This  dependence  exists  for 
the  inner  expansions  of  these  equations  as  well. 
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With  substitutions,  the  dfl/d£  equation  expansion  becomes 
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The  da/d£  Equation  Expansion.  It  has  been  noted  that 
the  da/d£  equation  is  a  function  of  the  dl/d£  equation  as 
well  as  a  function  of  a,  7,  and  I.  This  dependence 
exists  for  the  outer  expansions  of  these  equations  as  well 
From  Eq  (5.40) 


tana1 . dl1 


s  (1+e  £)tan7  tanl 
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With  substitution,  the  following  expression  is  obtained: 
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The  dq/df  Equation  Expansion.  The  outer  expansion  of 
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Inner  Expansion  Solutions 

e°  Terms .  The  complete  set  of  e°  term  inner  expansion 
differential  equations,  derived  above,  are  as  follows: 
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dIc 
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2u  Be 
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si 


1  L.cosa 


n70  CD  cos70 


CT  -  sinacosa 
L.r.-C. - - - o_ 


c°s70  sin70 


(5.45) 
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an  c,  , 


smasina. 


cos7  sin7  sinl 


CT  *  sinasina 

J±.  Be“^ - 5 _ 

CD  tanlosin7ocoslo 


CL  R 

7^ — *  Be  scosu 
CD 


(5.49) 


(5.51) 


(5.53) 


Solutions  to  this  set  of  differential  equations  are 
derived  in  Appendix  D  and  are  given  below.  Eq  (5.55)  is  the 
inner  zero  order  e  solution  to  the  dq/dh  equation  of  motion. 


'■'L  _  c 

=  Tr-'Be  ^cosa  +  K 


%  =  C 


_  y  2  2Cd7q 
Uo  3  %  ex^  C^cosa 


sina  sinl  =  sinK 

0  0  4 


cosaQ  =  cosK^  cos (K^  -  fi0 ) 


»IQ  =  cosK^  cos  tana  *  log  tan  (^  +  2^”)]  +  j 


(5.55) 


(5.56) 


(5.57) 


(5.58) 


(5.59) 


e 1  Terms .  The  complete  set  of  e1  term  inner  expansion 
differential  equations,  derived  above,  are  as  follows: 


2u  Be 
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cos70 

^1  .  2 
sin  7 


CT  CT  'i  cos  a 

+  ~  _L.cosa_ _ L.  _ 

*1  Cj.  COS7  .  2 

U  >o  D  sin7  cos  7 
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2u,Be“«-  -4—  *  ^.£282  .  0(e2) 

1  [sin70  CD  COS70J 


(5.46) 


Order  e  Inner  Solution  for  the  dq/dh  Equation  of  Motion. 


The  order  e  inner  solution  to  the  cosine  of  the  flight  path 
angle  is  given  by  the  substitution  of  Eqs  (5.55)  and  (5.60) 
into  Eq  (5 . 35) . 

q1  =  cJq  +  eqi  +  0(e2)  (5.35) 

CT 

q1  =  p^-Be  ^cosa  +  K  +  eK  (5.61) 

°D  1  7 

At  this  point,  the  inner  and  outer  solutions  to  the 

equations  of  motion  have  been  derived  to  order  e°  .  These 
are  the  solutions  to  the  atmospheric  entry  for  a  non¬ 
rotating  Earth.  In  addition,  the  inner  and  outer  solutions 
to  the  dq/dh  equation  of  motion  have  been  derived  to  order 
e.  These  latter  two  solutions  include  the  Coriolis 
acceleration  on  the  flight  vehicle  and  therefore  account  for 
the  significant  effects  of  a  rotating  Earth.  In  order  to 
create  composite  solutions  valid  for  both  the  inner  and 
outer  regions,  matching  is  performed. 

Matching  Zero  Order  e  Solutions 

One  of  the  fundamental  rules  of  the  technique  of  matched 
asymptotic  expansions  is  that  the  outer  expansion  of  the 
inner  expansion  solution  is  equal  to  the  inner  expansion  of 
the  outer  expansion  solution  (Nayfeh,  1981:277,278).  This 
condition  allows  matching  of  the  solutions,  reducing  the 
number  of  unknowns  in  the  solution  equations.  On  the 


.»  fct  l 


following  pages  the  inner  expansions  are  taken  of  the  e 


outer  expansion  solutions.  Next,  the  outer  expansions  are 


taken  of  the  e°  inner  expansion  solutions.  Finally,  the 


derived  expansions  of  expansions  are  matched  and  the  inner 


solution  constants  of  integration  are  expressed  in  terms  of 


the  outer  solution  constants  of  integration. 


Inner  Expansions  of  Order  £  Outer  Expansion  Solutions. 


From  Eq  (5.29),  the  one  term  outer  expansion  for  u  is 


(uQ  )  =  C/(l+h) 


2  2 

Noting  that  £  =  h/e  or  h  =  e  £  and  rewriting  this 


equation  in  terms  of  the  inner  variable,  f,  gii 


(u0)°  =  C/(l+e20 


Expanding  for  small  e  gives  the  inner  expansion  of  the 


outer  expansion  of  the  zero  order  e  term  for  u. 


K>°r- 


(5.62) 


From  Eq  (5.30),  the  one  term  outer  expansion  for  q  is 


(qo)°  =  C  /[2(1  +  h)  -  C  C2(l  +  h)2] 


Rewriting  this  equation  in  terms  of  the  inner  variable 


<V°  -  [v[2<» 


£)  -  C,C2  (1  +  e2{): 


Expanding  for  small  e  gives  the  inner  expansion  of  the  outer 
expansion  of  the  zero  order  e  term  for  q. 


K)T  -  C^/[2  *  -  C,°2  -  C,C2£<«  '  2C>C2£!<] 


it  A 


[(%)°Y  =  [c,/(2  -  CtC2)] 


1  /  2 


(5.63) 


From  Eq  (5.28),  the  one  term  outer  expansion  for  H  is 


Rewriting  this  equation  in  terms  of  the  inner  variable 
simply  gives 


=c, 


Expanding  for  small  e  gives  the  inner  expansion  of  the  outer 
expansion  of  the  zero  order  e  term  for  Q. 


K>°f  - 


(5.64) 


From  Eq  (5.27),  the  one  term  outer  expansion  for  I  is 


(I  )  =  C 

v  0  J  b 


Rewriting  this  equation  in  terms  of  the  inner  variable 
simply  gives 


Expanding  for  small  e  gives  the  inner  expansion  of  the  outer 


expansion  of  the  zero  order  e  term  for  I. 


[‘V0]1  -  °5 


(5.65) 


From  Eq  (5.31) ,  the  one  term  outer  expansion  for  a  is 


O  [1  -  C  /(1-h)  1 

'V  =  “  COB  (1-c’o'/’‘.  *  °3 

v  1  2  J 


Rewriting  this  equation  in  terms  of  the  inner  variable: 


o  f1  -  °,/(i+£  £)]  ,  c 

(“»)  =  ‘  “s  (1  -C»o  )>/•".  3 

v  1  2  J 

Expanding  for  small  e  gives  the  inner  expansion  of  the  outer 
expansion  of  the  zero  order  e  term  for  a. 


[  (ao  )  °]  =  ‘  COS 


1  -  C 

-1  _ 1  +  c 

.(1  -  C  2C  )'  /2.  3 

v  1  2  / 


(5.66) 


Outer  Expansions  of  the  e°  Inner  Expansion  Solutions 


From  Eq  (5.55),  the  one  term  inner  expansion  for  q  is 


(qQ)i  =  0^'Be  ^cosa  +  K] 


2  2  •  • 

Noting  that  £  =  h/e  or  h  =  e  £  and  rewriting  this 

equation  in  terms  of  the  outer  variable,  h,  gives 


cos?  +  K 


Expanding  for  small  e  gives  the  outer  expansion  of  the  inner 
expansion  of  the  zero  order  e  term  for  q. 


and  hence  )  *j  = 


(5.67) 


.,o 


cos  1  K 


(5.68) 


From  Eq  (5.56),  the  one  term  inner  expansion  for  u  is 

-^o 


(V1  =  kJ^)1; 


exp 


cosaC^/CpJ 

Rewritten  in  terms  of  the  outer  variable 

SCpCos'1  (c^)1 


(“o'1  =  K3  [<V']  exp 


C^cosa 


Expanding  for  small  e  gives  the  outer  expansion  of  the  inner 
expansion  solution  of  the  zero  order  e  term  for  u. 


-2cos  1 K 

.<U0M  =  K3K!  CXP 

l 

cosaC^/Cp 

(5.69) 


From  Eq  (5.59),  the  one  term  inner  expansion  for  I  is 


(I*)1  =  cos  1 


cosK  cos 
4 


tana  *  log  £ tan  ^ — jj  +  Kg 


Rewritten  in  terms  of  the  outer  variable 


(In ) 1  =  cos 


-1 


cosK  cos | tana 
4 


a)1- 


tana  *  log |tan ^  | — jj  +  Kg  j 


Expanding  for  small  e  gives  the  outer  expansion  of  the  inner 
expansion  of  the  zero  order  e  term  for  q. 


cos 


-1 


C  cos 
4 


cosK  cos  tana ' 


cos  1  K 


log [tan [*  ■*  2 - L]j  4  Ke] 


(5.70) 


From  Eq  (5.57),  the  one  term  inner  expansion  for  a  is 

sinK. 

(aj1  =  sin 


sin(IQ  )  1 


Rewritten  in  terms  of  the  outer  variable: 


sinK 


(a  l1  =  sin  1 

°  sin (I  ) XJ 


Expanding  for  small  e  gives  the  outer  expansion  of  the  inner 
expansion  of  the  zero  order  e  term  for  a. 


(5.71) 


r  •  i  ^ 

(«0)  J  =  sin 

sinK^ /sin[(IQ)  J 

wh 


ere 


|(IQ)*J  was  given  by  Eq  (5.65) 

From  Eq  (5.58),  the  one  term  inner  expansion  for  fi  is 
rcos (a  )  \ 

(00)1  =  -cos'1 [ " cosK "  ]  4  Kb 


Rewritten  in  terms  of  the  outer  variable 
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Vs 


(nA ) 1  =  -cos 


•  cos  (a  ) 1- 


Expanding  for  small  e  gives  the  outer  expansion  of  the  inner 


(5.72) 


expansion  of  the  zero  order  e  term  for  Q. 

.  ,  o  _  [  .  ,  o 

(no  )  1j  =  -cos  1  cos  (Oq)1]  /cosK4  +  K5 

f  il  ° 

where  [  (aQ )  I  was  given  by  Eq  (5.66) 


6°  Term  Matching.  As  previously  discussed,  the  outer 
expansion  of  the  inner  expansion  solution  is  equal  to  the 
inner  expansion  of  the  outer  expansion  solution . 


Mathematical ly , 


[wl]°  -  [«°f 

This  is  sometimes  expressed  (Lagerstrom,  1972:90)  as 

lim  ,  lim.  x  =  lim.  lim  .  x 
outL  in  J  ml  out  J 


(5.73) 


where  x=x  +  ex  +  e  x  + 

0  1  2 


The  expressions  derived  for  the  outer  expansion  of  the  inner 
expansion  solution  and  for  the  inner  expansion  of  the  outer 
expansion  solution  can  be  matched  for  q,  u,  a,  I,  and  Q. 

Matching  the  Speed  Ratio  Equations.  From  Eq  (5.73) 


K>f  -  k>°f 


Therefore,  from  Eqs  (5.62)  and  (5.69) 


C  =  K3K  exp 


-2cos  1  K 


[cosaCL/cDj 


Matching  the  Flight  Path  Angle  Equations 


Setting  Eq  (5.63)  equal  to  Eq  (5.57)  gives 


K  = 
1 


c, /(2  -  c  C2)] 


1  /  2 


Matching  the  Inclination  Angle  Equations 


[^>‘]°  -  [^.>°]X 


Equating  Eqs  (5.65)  and  (5.70)  gives 


C  =  cos 
5 


cosK  cos | tana 
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*  log 


cos  1  K 


V-/  O  XV 

l  +  — 2-4]  +  KeJ 


Matching  the  Argument  of  Latitude  at  Epoch 
Equations . 


K>r  ■  K>°r 


(5.74) 


(5.75) 


(5.76) 


Equating  Eqs  (5.66)  and  (5.71)  gives 


(1  -  c  2c2)l/2J  +  °3 


iin  sinK^ /sin[(Io)  J 


and  since  [(I^1]  =  [  CI0  ^  °]  =  C5  th< 


1  -  C 


1  ,  o  •  -1  4 

- 7 —  +  C  =  sm  — ; pr~ 
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sinK. 


(5.77) 


Matching  the  Longitude  of  the  Ascending  Node 


Equations 


K)1]’  ■  K>°f 


Setting  Eq  (5.64)  equal  to  Eq  (5.72) 


_  r  J.  . ,  o 

C  =  -cos  1  cos  (a  ) 1  /cosKi  +  K 
A  L  0  J  4  5 


r  l  -  c 


C  =  -cos  cos  -  COS  - J- 

4  [  [  1(1  -  c/cy’/2 


+  C  / cosK  +  K 

3  '  4  S 


(5.78) 


The  inner  expansion  solutions,  valid  in  the  domain  of  h 


near  the  boundary  layer  where  aerodynamic  forces  are 


dominate,  have  been  matched  with  the  outer  expansion 


solutions,  valid  in  the  domain  of  h  away  from  the  boundary 


layer.  The  results  of  the  matchings  are  used  to  reduce  the 


number  of  unknowns  in  the  solution  equations.  The  constants 


of  the  outer  expansions  are  determined  first,  often  from 


inj™r»wyxv*-»>  ,-*  , 


•>  •jr’yry  v  ’.y  v  v/.’/v.v.v.v;.’.'v.v.v.v.  v.t.v.-v 


initial  conditions,  in  a  planetary  entry  problem.  Hence 
is  desirable  to  express  the  constants  of  the  inner 
expansions,  K,  in  terms  of  the  constants  of  the  outer 
expansions,  C.  From  Eq  (5.75) 

K,  -  [V<2  -  °,C2>],/!  <5 

From  Eqs  (5.74)  and  (5.75) 


K  =  (2  -  C  C 

3  v  12 


-2cos  1  [C5 /(2  -  C  C2)] 
)exp  cosa • CL/CD 


From  Eq  (5.77) 


K  =  sin  1  sinC  *sin  -  cos  5 


1  -  C, 

- - - 7-  +  C 

(1  -  C  2C  )l/2.  3 

12 


From  Eq  (5.78) 


.  1  -  C 

V  n  .  1  1 

K  =  C  +  cos  COS  -  COS  - T~ 

54  [  l  [(1  -  c  2c2)l/z 


+  C  /cosK 
3  4 


From  Eq  (5.76) 


cosC 
“I  5 

K  =  cos  — — jj— 

6  cosK 

4 


cos  1  K  , 

• 

*  .  1 

tan 

[4  2  J 

With  Eqs  (5.79)  -  (5.83)  the  inner  solutions  can  be 


written  in  terms  of  the  constants,  C. 


Zero  Order  e  Solutions  to  the  Equations  of  Motion 

A'  solution  good  for  the  entire  range  of  h  is  desired  for 
the  convenience  of  not  having  to  apply  the  two  sets  of 
solutions  to  one  entry  problem.  The  composite  solution, 
valid  everywhere  in  the  domain,  combines  the  inner  and  outer 
solutions  that  were  valid  only  in  certain  overlapping 
regions  of  the  non-dimensional  altitude  domain.  This 
composite  solution  is  given  in  the  following  equation  by  the 
method  of  matched  asymptotic  expansions. 

(y)°]  =  y°  +  y1  -  [(y)1]  (5.84) 

The  composite  solution  is  therefore  given  by  the  sum  of 
the  inner  and  outer  solutions  and  the  difference  between 
this  result  and  any  expression  common  to  both  the  inner  and 
outer  solutions.  In  addition,  the  inner  expansion  of  the 
composite  solution  is  equal  to  the  inner  solution,  and  the 
outer  expansion  of  the  composite  solution  is  equal  to  the 
outer  solution.  Mathematically, 


and 


(5.85) 


Speed  Ratio  Composite  Solution . 


c  /  \  °  /  .  i 

u  =  (u)  +  (u) 


(u)°  +  (u)1 
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•T-  *  V 


Combining  Eqs  (5.29),  (5.56),  and  (5.62)  gives  the  composite 
solution  for  the  speed  ratio. 


=  °,MhI  -  >]  *  K3[|-Be'h/£  ♦  k,]2- 

•exp  -2cos  *  Be  '  •  cosa  +  K  /^cos  a 

lCD  1 J  °D  J 


Flight  Path  Angle  Composite  Solution, 


(5.86) 


q°  -  (q)°  +  (q)1  ~  [  Cq)°]  1  =  (q)°  +  (q) 1  -  [  (q) £]  ° 

Combining  Eqs  (5.30),  (5.55),  and  (5.67)  gives  a  composite 

solution  for  the  cosine  of  the  flight  path  angle. 


qc  =  Cj  / [2 (1+h)  -  C^U+h)2] 


CL  „  -h/e2 
+  p— -Be  '  1 


(5.87) 


Inclination  Angle  Composite  Solution. 


Ic  =  (I)°  +  (I) *  -  [(I)0]1  =  (I)°  +  (I)1  -  [  CD  i]° 

Combining  Eqs  (5.27),  (5.59),  and  (5.70)  gives  a  composite 

solution  for  the  orbital  inclination  angle. 


C  “1 

I  =  cos  cosK  cos  tana ' 


(5.88) 


,  .  k  1  -1  Ln  -h/e2  v 

log  tan  -r  +  ^cos  p— Be  cosa  +  K 

4  z  '■'n  1 


Combining  Eqs  (5.31),  (5.59),  and  (5.66)  gives  a  composite 

solution  for  the  longitude  of  the  ascending  node. 
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(1  -  C  2C  )'  /2 
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(5.89) 


Argument  of  Latitude  at  Epoch  Composite  Solution . 

nc  =  (n)°  +  (n)1  -  [(n)0]1  =  (n)°  +  (n)1  -  [ (n) i] 

Combining  Eqs  (5.28),  (5.58),  and  (5.72)  gives  a  composite 
solution  for  the  argument  of  latitude  at  epoch. 


fi  =  -cos 


rcos(o1)1 

[  cosK  ]  +  KS 
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nc  =  K,  -  cos  5 


cos 


sin  1  j^sinK^ /sin  (IC)  j 


/cosK 


(5.90) 


For  Eqs  (5.86)  -  (5.90),  the  constants  are  given  by 
Eqs  (5.79)  -  (5.83) .  Hence,  solutions  to  the  five  equations 


of  motion  have  been  developed  as  functions  of  , 

C  ,  and  C  . 

4  5 


Matching  £  Order  Solutions  to  the  dq/dh  Equation  of  Motion 
An  inner  expansion  of  the  outer  expansion  solution  is 
undertaken  as  well  as  an  outer  expansion  of  the  inner 
expansion  solution.  Matching  is  performed  by  equating  these 
expansions  of  expansion  solutions. 

Inner  Expansion  of  the  Outer  Expansion  Solution.  The 
outer  expansion  solution  for  the  dq/dh  equation  of  motion 
was  derived  to  order  e  and  given  by  Eq  (5.34) . 
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2 (1+h) 
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-  c2  (1+h)2 


-1  /  2 


(C  )2 
v  1  6  J 


C20e*»[-C,6h] 


2  2 

Noting  that  £  =  h/e  and  h  =  e  £  and  rewriting  the 
above  equation  in  terms  of  the  inner  variable,  £,  gives 

I*  A 
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2(1  +  e!0  -  (1  ■.  €25)2] 


-  e 


13  1  4 


1  6 


(C  ) 
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-C,6£ 


<] 


(5.81) 


Expanding  for  small  e  gives  the  inner  expansion  of  the  outer 
expansion  solution  to  order  e . 


[(q)0]1  =  [C./2  -  C  C,)] 


C1  3  _  Ci  4  _  c 

Cl 6  (c6)2  20 


(5.92) 


Outer  Expansion  of  the  Inner  Expansion  Solution .  The 
inner  expansion  solution  for  the  dq/dh  equation  of  motion 
was  derived  to  order  e  and  given  by  Eq  (5.61) . 


(q)  *  =  pr^-Be  ^cos a  +  K  +  eK 
CD  1  7 


Rewriting  this  equation  in  terms  of  the  outer  variable, 
h,  and  expanding  for  small  e  gives  the  outer  expansion  of 
the  inner  expansion  solution  for  q. 


[(q)1]0  =  K  -  eK? 


(5.93) 


Matching .  The  outer  expansion  of  the  inner  expansion 
solution  is  equated  to  the  inner  expansion  of  the  outer 
solution  expansion. 


[(q)f  -  [(q)°f 


11/2  fC,3  CH 

K  +  eK  =  C  /2  -  C  C  )  -  e  - x—  -  C 
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Subtracting  Eq  (5.79)  from  this  expression  and  dividing  by  e 


gives 


The  composite  solution  to  the  dq/dh  equation  can  be 


easily  constructed: 

qc  =  (q)°  +  (q)1  -  (q)1]  =  (q)°  +  (q)1  -  [(q)°] 
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(5.95) 


where  C  ,  C  ,  and  C  „  are  given  by  the  following  relations 
1314'  16  &  & 


derived  in  Appendix  C: 


C,  3  =  2cos(Cb) 
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16  ”  c,  (c  i),/2  C.CU 


-  1 


(5.98) 


where 


and  where  C  and  C  are  constants  of  integration. 

6  2  0 


Hence,  a  solution  for  q  has  been  found  as  a  function  of  two 


constants  of  integration,  C  and  C  „ ,  and  the  constants  C  , 

&  ’6  20’  13’ 


'  ,  and  C  ,  which  are  dependent  upon  C  and  C  . 

14  16  r  1  2 


Summary 


For  convenience,  a  summation  of  the  derivation  results 
using  matched  asymptotic  expansions  are  presented. 


Solutions  to  the  Non-Rotating  Earth  Equations  of  Motion , 


Treatment  of  atmospheric  entry  as  a  boundary  layer 
problem  allowed  for  the  application  of  the  method  of  matched 
asymptotic  expansions.  In  this  section,  the  rotating  Earth 
equations  of  motion  were  asymptotically  expanded  using  a 
small  parameter  that  is  a  function  of  the  planet’s  inverse 
atmospheric  scale  height,  radius,  and  rotation  rate.  Zero 
order  solutions  to  these  equations  were  derived.  These 
solutions  were  of  two  types,  inner  solutions,  valid  close  to 
the  Earth,  and  outer  solutions,  valid  far  from  the  Earth. 

The  domains  of  these  two  sets  of  solutions  overlapped. 

Hence,  direct  matchings  of  the  expansions  of  the  expansion 
solutions  were  accomplished.  This  produced  equations  that 
related  the  inner  solution  constants,  K. ,  to  the  outer 

solution  constants,  .  Composite  solutions,  valid  for  the 

entire  atmospheric  entry  domain  of  non-dimensional  altitude, 
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were  constructed  using  these  relations.  These  composite 
solutions  are  solutions  to  the  non-rotating  Earth  equations 
of  motion  for  three-dimensional,  lifting  atmospheric  entry. 


The  outer  solutions  were  given  by 


I  =  C 
0  s 


=  C„ 
0  4 


u  (1+h)  =  C 
o  v  '  1 


^  -  V1+h)2 


ri  -  c  / (i+h)  ■ 

a  —  -  cos  - 7—  +  C 

0  .  (1  -  C  2C  )’  /2.  3 

v  1  2  ' 


(5.27) 


(5.28) 


(5.29) 


(5.30) 


uo  =  1  +  cos  (a0  -  03)[l  -  c/cj 


2n  1  1  /2 


(5.31) 


The  inner  solutions  were  given  by 


W  _  c 

a  =  pr~'Be  ^  coso  +  K 
^  °D  1 


K  2  CD^0 

u  =  K  a  exp  ^ - - 

o  3  \)  r  C^cosa 


sina  sinl  =  sinK 

o  o  A 


cosaQ  =  cosK^ cos (K&  -  ) 


(5.55) 


(5.56) 


(5.57) 


(5.58) 


/y 

cosIQ  =  cosK^  cos  tana ‘log  tan  (^  +  2^]  +  ^6 


(5.59) 


Matching  the  expansions  of  these  solution  sets  produced  the 

following  relationships  between  the  outer  solution 

constants,  C. ,  and  the  inner  solution  constants,  K. . 

’  1  ’  ’  1 
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K  =  sin 
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-1 


sinC  ’sin 
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1  -  c 
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5 

3 
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/ cosK^ 
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cosC 

0 

cosK 

A 

(5. 

The  equations  given  above  are  for  the  three-dimensional 
atmospheric  entry  of  a  non-rotating  Earth.  These  solutions 
were  derived  from  zero  order  terms  in  the  expansions  of  the 
rotating  Earth  equations  of  motion. 

The  composite  solutions  for  the  three-dimensional 
atmospheric  entry  of  a  non-rotating  Earth  were  given  by 


tana • log 


r  cos  K  , 

tan 

*  4  .  ’ 

L4  2  J 

K  =  cos 
6 


lie 


u  =  SlxiThy  -  JJ  +  K3lcfBe  '  -cos*  +  K,J 


exp  -2cos  1  IT^-Be  h/e2‘ 
L  l) 


cosa  +  K  /jsr-cos  a 
:J  CD 
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qc  =  Cj / [2 (i+h)  -  c  c2  (i+h)2] 


CD 


(5.87 


1°  =  cos  1  cosK  cos  tana 
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nc  =  K  -  cos  1 
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cos  sin  1  fsinK  /sin (IC) 1  /cosK  (5.90 
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Solution  to  the  dq/dh  Rotating  Earth  Equation  of  Motion 


To  accurately  model  Earth  atmospheric  entry,  a  solution 


to  the  dq/dh  equation  that  included  rotating  Earth  effects 
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was  required.  This  solution  was  derived  in  the  form  of 


order  e  inner  and  outer  solutions  to  the  asymptotically 


expanded  dq/dh  equation  of  motion  for  a  rotating  Earth. 


Matching  was  performed  between  these  solutions  and  a 


composite  solution  to  the  dq/dh  equation  of  motion  for  the 


three-dimensional  atmospheric  entry  of  a  rotating  Earth  was 


constructed . 


The  outer  solution  to  the  rotating  Earth  dq/dh  equation 


of  motion  was  given  by 


o 

q  = 


-  c2(i+h) 


1  “I  / 2 


rc  c  c 

1  3  _  1  4  _  ,  14 

'  <C16)2  '  C1  6 


C  exp  -C  h] 
2  0  r  1  6  J 


(5.34) 


The  inner  solution  to  the  rotating  Earth  dq/dh  equation  of 


motion  was  given  by 


q*  =  ^r^'Be  ^cos o  +  K  +  eK 
°D  1  1 


(5.61) 


The  composite  solution  to  the  dq/dh  equation  of  motion  is  of 


order  e  and  hence  includes  Earth  rotational  terms.  The 


composite  solution  to  the  dq/dh  equation  of  motion  for  the 


three-dimensional  atmospheric  entry  of  a  rotating  Earth  was 


given  by 


qc  =  Cj / [2 (i+h)  -  cc2  (1+h)2] 
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where  C  „ ,  C  ,  and  C  were  given  by  the  following 
1314  16  0  & 

relations  derived  in  Appendix  C: 


C1  3  =  2cOS  (CS  >  [c^  -  'J  '  IC7C77 


(Cj  )2  (0t  J  )3  /2 


3C  C  C  3 

C  =  - -  6  12  ,  +  2cos  (C  )  1  -  3-^-  •  — - - 
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where 


CJ  and  c,2  =  [t-cJ 


and  where  and  q  are  constants  of  integration 
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Traiectorv  States  of  Validity  for  the  Non-Rotatin 


Earth  Equations  of  Motion 


Introduction 

In  Section  IV,  the  rotating  planet  terms  in  each  of  the 
five  equations  of  motion  for  Earth  atmospheric  entry  were 
examined.  The  dq/dh  equation  of  motion  for  a  non-rotating 
Earth  was  found  to  be  never  valid  for  the  investigated 
ranges  of  inclination  angle  and  speed.  Because  of  this,  a 
solution  was  developed  for  the  dq/dh  equation  of  motion  for 
the  three-dimensional  lifting  atmospheric  entry  of  a 
rotating,  spherical  Earth.  This  was  accomplished  in  Section 
V  by  treating  atmospheric  entry  as  a  boundary  layer  problem 
and  applying  the  method  of  matched  asymptotic  expansions. 

In  Section  IV,  trajectory  states  were  found  to  exist 
where  some  of  the  non-rotating  Earth  equations  of  motion  are 
valid  for  a  rotating  Earth.  The  d(l/dh,  da/dh,  and  dl/dh 
non-rotating  equations  of  motion  were  all  found  to  be  valid 
for  the  same  entry  trajectory  states.  Other,  independent 
trajectory  states  were  found  to  exist  where  the  du/dh  non¬ 
rotating  equation  of  motion  is  valid  for  a  rotating  Earth. 

In  this  section,  these  trajectory  states  are  examined  in 
more  detail  for  the  du/dh  and  dl/dh  equations  of  motion. 
Plots  of  the  Rotate  term  solutions  to  these  equations  are 
generated  for  a  large  range  of  realistic  values  of  u,  h,  a, 
7,  and  I.  Trends  in  these  solutions  are  discussed  and  the 


non-existence  of  overlapping  trajectory  state  solutions 
between  the  du/dh  and  dl/dh  Rotate  equations  is  verified, 
few  methods  to  estimate  solutions  to  the  du/dh  and  dl/dh 
Rotate  equations  are  also  presented. 

Solutions  to  the  du/dh  Rotate  Equation.  In  Section  IV 
it  was  found  that  the  rotating  planet  terms  in  the  du/dh 
equation  of  motion  are 


I  Rotate 


l1  /2 

=  -  4  e2u(l+h)|  cosl 


-  2e2 (1+h) 


.  2  T  . 

2  cosasin  lsinq 
tan^ 


(6.1) 


Setting  this  equation  equal  to  zero  and  solving  for  flight 
path  angle  gives  solutions  for  trajectory  states  where  the 
rotating  Earth  terms  in  the  du/dh  equation  of  motion  are 


7  =  tan 

Rotate  =  0 


2 

-l  -Ccosqsinasin  I 
2cosI 


(6.2) 


where  C  =  1 

u  2  J 


(6.3) 


Solutions  to  the  dfi/dh.  da/dh,  and  dl/dh 


Rotate  Equations .  In  Section  IV,  the  solutions  for  the 
Rotate  term  expressions  were  found  to  be  identical  for  the 
da/dh,  dfl/dh,  and  dl/dh  equations  of  motion.  This  result 
simplifies  the  search  for  trajectory  states  where  the  non¬ 
rotating  Earth  equations  of  motion  are  valid  for  a  rotating 


Earth.  Instead  of  detailed  examination  of  the  solutions  to 


all  three  of  the  da/dh,  dn/dh,  and  dl/dh  Rotate  term 
expressions,  examination  of  only  one  of  them  is  required. 
In  this  section  the  solutions  to  the  dl/dh  Rotate  term 
expression  are  examined.  From  Eq  (4.13) 


I  Rotate 


=  2  fe  iim 
2  u 


cosQsinl 

tan^ 


(cosatan^  -  sina) 


,,  ,  .2  cosasinacosisinl 
-  e  (1+h)  * - x - 

2  utan^ 


(6.4 


Setting  this  equation  equal  to  zero  and  solving  for  flight 
path  angle  gives  solutions  for  trajectory  states  where  the 
rotating  Earth  terms  in  the  dl/dh  equation  of  motion  are 


.  -l  sina  +  Ccoslsina 

7  =  tan  - o - 

Rotate  =  0  Zcosa 


(6.5; 


Non-Existence  of  Overlapping  Solutions  to  the  du/dh  and 


dl/dh  Rotate  Equations 

To  verify  the  non-existence  of  overlapping  solutions  to 
the  du/dh  Rotate  and  dl/dh  Rotate  equations,  Eqs  (6.2)  and 
(6.5)  are  equated. 


-1  -Ccosasinasin  I 


-1  sina  +  Ccoslsina 


A  trivial  solution  to  this  expression  is  given  by  sina  =  0  . 
The  non-trivial  solutions  to  Eq  (6.6)  are  found  by  dividing 
by  sina  and  then  solving  for  cosa. 


f-cosl  2  X1  1  /  •  T 

cosa  =  I  —  -  cos  II  /sinl 

„h„e  c  = 


(6.7) 


For  real  solutions  to  exist  for  Eq  (6.7),  the  following 
condition  must  be  true: 


cosl 

C 


+ 


cos2 I  <  0 


(6.8) 


Since  values  of  non-dimensional  altitude  and  speed  ratio 
are  never  negative,  C  always  has  a  positive  value.  The  term 
cosl  is  never  negative  for  prograde  entry  trajectories. 
Hence,  the  condition  given  by  Eq  (6.8)  is  never  satisfied 
and  overlapping  solutions  to  the  du/dh  and  dl/dh  Rotate 
equations  do  not  exist. 


Graphical  Trajectory  State  Examination 


combinations  of  u,  h,  a,  and  I  into  the  Rotate  solution 
equations,  given  by  Eqs  (6.2)  and  (6.5),  for  the  du/dh  and 
dl/dh  equations  of  motion.  These  solutions  give  trajectory 
states  for  which  the  non-rotating  equations  of  motion  are 
perfectly  valid  for  rotating  Earth  entry  analysis.  This 
validity  is  dependent  upon  the  physical  characteristics  of 
the  planet  under  consideration,  namely  e?  ,  and  the  type  of 

trajectory  being  analyzed;  validity  is  not  directly 
dependent  on  any  of  the  vehicle’s  physical  characteristics. 
Determination  of  the  validity  of  the  non-rotating  equations 
of  motion  for  planets  other  than  Earth  could  be  undertaken 
in  a  similar  manner  as  the  analysis  given  in  this  section. 

A  computer  program  was  developed  to  solve  Eqs  (6.2)  and 
(6.5)  for  a  large  range  of  u,  h,  a,  and  I.  Plots  of 
families  of  curves  of  7  versus  I  are  presented  for  various 
ranges  of  a  and  for  realistic  combinations  of  u  and  h. 
Appendix  E  presents  three  example  Earth  atmospheric  entry 
trajectories  to  give  an  indication  of  what  the  dimensionless 
variables  u  and  h  are  in  relation  to  more  conventional  entry 
trajectory  parameters  such  as  Mach,  altitude,  and  velocity. 
Most  current  and  planned  lifting  entry  vehicles  have 
trajectories  which  fall  within  the  ranges  of  u  and  h  given 
by  these  examples.  The  realistic  combinations  of  u  and  h 
discussed  in  the  following  pages  are  derived  from  liberal 
estimation  of  their  ranges  in  these  example  trajectories. 
Values  of  u  selected  for  investigation  ranged  from  circular 


orbital  velocities  to  low  supersonic  speeds  where  terminal 
maneuvers,  such  as  landing  approaches,  are  usually 
initiated.  Realistic  values  of  non-dimensional  altitude 
corresponding  to  these  values  of  speed  ratio  were  used  in 
this  analysis.  Values  of  argument  of  latitude  at  epoch  and 
orbital  inclination  angle  had  ranges  of  0.0  <_£»<_  360 
degrees  and  0 . 5  I  £  75  degrees,  where  most  atmospheric 
entry  occurs.  Because  of  the  large  number  of  plots  needed 
to  display  data  trends,  some  of  the  plots  which  are 
discussed  in  this  section  are  contained  in  Appendix  F.  In 
all  the  plots  referred  to  in  this  section,  a,  the  argument 
of  latitude  at  epoch,  is  referred  to  as  "Alpha".  All  values 
of  7,  a,  and  I  given  in  these  plots  are  in  units  of  degrees. 

Trajectory  States  of  Validity  for  the  du/dh  Non-Rotating 

Earth  Equation  of  Motion.  Trajectory  states  are 
presented  below  where  the  non-rotating  Earth  du/dh  equation 
of  motion  is  valid.  Trends  in  these  solutions  to  Eq  (6.2) 
are  discussed. 

IZR  ical  Solution  Observations.  Figure  13  is  a 
typical  plot  of  solutions  to  Eq  (6.2)  of  flight  path  angle 
versus  orbital  inclination  angle.  For  constant,  non-zero 
values  of  argument  of  latitude  at  epoch,  7  decreases  as  I 
increases.  Values  of  flight  path  angle  also  decrease  in 


this  plot  as  a  increases.  In  addition,  as  a  increases  for  a 
constant  I,  the  change  in  7  decreases,  especially  for  large 
inclination  angles. 


13. 


Effect  of  h  Variation.  Holding  all  other  variables 


constant,  the  effect  of  change  in  non-dimensional  altitude 
on  these  solutions  is  negligible.  Figures  FI  and  F2  in 
Appendix  F  present  results  for  values  of  h  higher  and  lower 
than  in  Figure  13  for  the  same  values  of  u  and  a. 

Comparison  of  Figures  FI  and  F2  indicates  that  only 
extremely  small  changes  in  flight  path  angle  occur  for  this 
large  change  in  h.  This  ineffectiveness  of  h  variation  on  7 
can  also  be  intuitively  seen  from  examination  of  Eq  (6.3). 
Note  that  for  all  realistic  values  for  non-dimensional 
altitude,  1.0  is  much  larger  than  h  in  the  expression 

a  *  h)3/2 

Effect  of  u  Variation.  Figures  F3  and  F4  present 
solution  results  for  values  of  speed  ratio  higher  and  lower 
than  in  Figure  13  for  the  same  range  of  values  of  argument 
of  latitude  at  epoch.  Comparison  of  Figures  13,  and  Figures 
FI  -  F4  demonstrates  that  large  changes  in  flight  path  angle 
occur  for  this  large  range  of  u.  As  u  decreases,  the 
magnitude  of  7  increases  for  curves  of  constant  a 

Effect  of  a  Variation.  The  trend  of  the  effect  of 
variation  in  argument  of  latitude  at  epoch  can  be  seen  by 
review  of  F.gures  13  and  Figures  FI  -  F4.  For  the  range  of 
values  of  0  in  these  plots,  7  and  the  change  in  7  decrease 
as  a  increases  along  lines  of  constant  inclination  angle. 
These  trends  are  more  prominent  for  large  values  of  I.  The 
a  terms  in  Eq  (6.2)  are 
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-  cos(a)sin(a) 


Examination  of  this  term  and  the  properties  of  the  inverse 


tangent  function  in  Eq  (6.2)  indicates  the  solutions  given 


in  Figure  13  will  be  symmetric  about  a  =  0  and  180  degrees 


and  these  symmetric  solutions  will  repeat  about 


±90  -  a  (degrees)  for  ±  a 


Figures  13,  F5,  F6,  and  F7  demonstrate  these 


observations . 


Trajectory  States  in  Three  Dimensions.  Figures  14  and 


15  present  three-dimensional  plots  of  the  trajectory  states 


of  validity  for  the  du/dh  non-rotating  Earth  equation  of 


motion.  Solutions  to  Eq  (6.2)  are  plotted  here  for  flight 


path  angle  versus  orbital  inclination  angle  and  argument  of 


latitude  at  epoch.  The  solution  surfaces  in  these  figures 


are  constructed  of  lines  of  constant  a  and  lines  of  constant 


flight  path  angle.  These  contours  of  constant  7  are  drawn 


for  changes  of  7  of  0.165  degrees  in  Figure  14  and  changes 


of  7  of  1.04  degrees  in  Figure  15.  In  both  of  these  plots, 


values  of  flight  path  angle  are  near  zero  for  small 


inclination  angles  and  for  very  small  values  of  argument  of 


latitude  at  epoch.  Although  Figures  14  and  16  present 


solutions  to  Eq  (6.2)  for  vastly  different  values  of  speed 


ratio,  the  two  plots  look  almost  identical;  only  the  scale 


of  flight  path  angle  values  differs. 
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'  States  of  Validity  for  the  dl/dh  Non-Rotatin 


Earth  Equation  of  Motion.  Trajectory  states  are  given  below 
where  the  non-rotating  Earth  dl/dh  equation  of  motion  is 
valid.  Trends  in  these  solutions  to  Eq  (6.5)  are  discussed. 

Typical  Solution  Observations.  Figure  16  is  a 
typical  plot  of  solutions  of  flight  path  angle  versus 
orbital  inclination  angle.  For  constant  non-zero  values  of 
argument  of  latitude  at  epoch,  7  increases  as  I  increases. 
Values  of  flight  path  angle  decrease  in  this  plot  as  values 
of  a  increase . 


Effect  of  h  Variation.  Holding  all  other  variables 
constant,  the  effect  of  change  in  non-dimensional  altitude 
on  these  solutions  is  negligable.  Figures  F8  and  F9  present 
solution  results  for  values  of  h  higher  and  lower  than  in 
Figure  16,  and  for  the  same  values  of  u  and  a.  Comparison 
of  Figures  F8  and  F9  indicates  that  only  extremely  small 
changes  in  7  occur  for  this  large  change  in  h.  As  in  the 
case  for  the  du/dh  Rotate  solutions,  this  ineffectiveness  of 
h  variation  on  7  can  also  be  intuitively  seen  from 
examination  of  the  expression  for  C  in  Eq  (6.5) . 

Effect  of  u  Variation.  Figures  F10  and  Fll  present 
solutions  for  values  of  speed  ratio  higher  and  lower  than  in 
Figure  16  for  the  same  range  of  values  of  argument  of 
latitude  at  epoch.  Comparison  of  Figures  16,  F10,  and  Fll 
indicates  that  large  changes  in  7  occur  for  this  large  range 
of  u.  As  u  decreases  (along  with  the  ineffectual  non- 
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dimensional  altitude),  values  of  flight  path  angle  increase 


for  curves  of  constant  a. 


Effect  of  a  Variation.  The  trend  in  the  effect  of 


variations  of  argument  of  latitude  at  epoch  can  be  seen  by 


review  of  Figures  16  and  Figures  F8  -  Fll.  For  the  range  of 


a  in  these  plots,  7  decreases  as  a  increases  for  constant  I. 


This  change  in  flight  path  angle  is  more  prominent  for  small 


values  of  speed  ratio.  Change  in  the  range  of  a  can  also  be 


easily  seen  analytically.  The  a  terms  in  Eq  (6.5)  can  be 


written  simply  as  tan (a)  .  Examination  of  this  term  and 


the  properties  of  the  inverse  tangent  function  in  Eq  (6.5) 


indicates  the  solutions  given  in  Figure  16  will  be  symmetric 


about  a  =  0  and  180  degrees,  and  that  these  symmetric 


solutions  will  repeat  about 


q  *  180  (degrees)  for  *  a 


Figures  16,  F12,  F13,  and  F14  demonstrate  these 


observations . 


Trajectory  States  in  Three  Dimensions.  Figures  17,  18, 


and  19  present  three-dimensional  plots  of  the  trajectory 


states  of  validity  for  the  dl/dh  non-rotating  Earth  equation 


of  motion.  Solutions  to  Eq  (6.5)  are  plotted  here  for 


flight  path  angle  versus  orbital  inclination  angle  and  for 


argument  of  latitude  at  epoch.  The  solution  surfaces  in 


these  figures  are  constructed  of  lines  of  constant  flight 


path  angle  and  lines  of  constant  argument  of  latitude  at 


epoch.  Figures  17,  18,  and  19  display  plots  generated  for 


^  * 
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Figure  18.  Surface  of  Trajectory  State  Solutions  of  Validity  for  the 
Rotating  Earth,  dl/dh  Equation  of  Motion  (u  =  .300  ,  h  = 


wvtot; 


successively  decreasing  values  of  speed  ratio.  Although 
these  values  of  u  differ  greatly,  the  solution  surfaces  are 
very  similar  in  appearance. 


Solution  Estimation  Methods 

Some  effort  was  spent  in  developing  the  computer  program 
to  solve  Eqs  (6.2)  and  (6.5)  for  large  ranges  of  values  of 
speed  ratio,  non-dimensional  altitude,  orbital  inclination 
angle,  flight  path  angle,  and  argument  of  latitude  at  epoch. 
To  rapidly  assess  the  validity  of  the  non-rotating  Earth 
atmospheric  entry  equations  of  motion  for  a  particular 
trajectory  state,  there  is  an  advantage  to  having  a  quick 
estimation  algorithm.  On  the  following  pages,  methods  are 
presented  to  estimate  some  solutions  to  the  du/dh  and  dl/dh 
Rotate  expressions  given  by  Eqs  (6.2)  and  (6.5), 
respectively . 

Estimates  to  the  du/dh  Rotate  Solutions.  Figure  20 
presents  a  plot  of  trajectory  state  solutions  to  Eq  (6.2), 
the  du/dh  Rotate  equation.  The  solution  surface  in  this 
figure  is  onstructed  of  lines  of  constant  flight  path  angle 
and  lines  of  constant  argument  of  latitude  at  epoch.  Two 
separate  estimations  of  solutions  to  Eq  (6.2)  can  be  made 
with  plots  such  as  Figure  20.  The  first  is  the 
approximation  that  values  of  flight  path  angle  in  the 
triangular  plane  (ABC)  of  trajectory  states  have  some 
constant,  near-zero  value.  In  general,  this  is  a  good 
approximation  for  7,  especially  for  high  values  of  u.  The 


ABC  plane  of  "constant"  7  is  bounded  by 


0.0  <_  a  <_  45.0  degrees 
0.5  <  I  75.0  degrees 

Line  AB :  a  =  -(0.60403)1  -  (45.302)  degrees 

Therefore,  the  conditions  for  assuming  flight  path  angle  has 
some  near-zero  value  are  given  by 

( . 02207) a  +  (.01333)1  <  1.00  (6.9) 

and  by  the  ranges  of  inclination  angle  and  argument  of 
latitude  at  epoch  given  above. 

The  maximum  error  in  assuming  7  =  0.0  degrees  in  the 
ABC  plane  in  Figure  20  is  approximately  equal  to  the  value 
of  7  corresponding  to  the  second  flight  path  angle  contour 
from  the  top  of  the  plot.  This  maximum  error  in  flight  path 
angle  is  about  7  =  -0.56  degrees  for  this  plot. 

Application  of  this  estimate  to  the  corresponding  ABC  planes 
in  Figures  14  and  15  give  maximum  errors  in  flight  path 
angle  of  about  0.30  and  1.5  degrees,  respectively. 

An  estimate  of  a  non-zero  value  for  flight  path  angle  in 
the  ABC  plane  can  reduce  these  maximum  errors.  Assuming  the 
value  of  7  in  this  plane  is  equal  to  half  the  7  value  of  the 
second  constant  7  contour  in  Figure  20  produces  an  estimated 
flight  path  angle  of  -0.28  degrees  with  maximum  error  of 
0.28  degrees.  Likewise,  estimates  of  7  =  -0.15  degrees 
and  7  =  -0.75  degrees  in  the  ABC  planes  in  Figures  14  and 
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15  give  maximum  errors  in  flight  path  angle  of  about  0.15 
degrees  and  0.75  degrees,  respectively. 

Another  approximation  to  the  solution  of  Eq  (6.2)  is 
given  by  the  simple  assumption  that  flight  path  angle  has 
near-zero  values  for  small  values  of  inclination  angle. 
Maximum  error  in  an  estimated  flight  path  angle  is  dependent 
upon  the  chosen  value  of  inclination  angle  for  this 
approximation.  For  example,  assuming  7  =  0.0  degrees  for 
I  15.4  degrees  invokes  a  maximum  error  in  flight  path 
angle  of  about  0.28  degrees  in  Figure  20.  Assuming  7  =  0.0 
degrees  for  I  29.2  degrees  produces  a  maximum  error  in 
flight  path  angle  of  about  0.56  degrees  in  Figure  20. 

An  estimate  of  a  non-zero  value  for  flight  path  angle  in 
the  ABC  plane  can  reduce  these  maximum  errors  for  the  same 
range  of  inclination  angle.  Assuming  7  =  -0.14  degrees 
for  I  <  15.4  degrees  produces  a  maximum  error  in  7  of 
about  0.14  degrees  in  Figure  20.  Estimation  of  a  flight 
path  angle  value  of  -0.28  degrees  for  I  29.2  degrees 
gives  a  maximum  error  of  0.28  degrees. 

The  accuracy  of  both  of  these  estimations  for  solutions 
to  the  du/dh  Rotate  equation  is  heavily  dependent  on  the 
value  of  speed  ratio  for  a  given  trajectory  state  and  the 
estimate  of  flight  path  angle. 

Planar  Estimation  to  dl/dh  Rotate  Solutions.  A  planar 
estimation  method  was  developed  to  obtain  rough  estimates  of 
solutions  to  Eq  (6.5).  For  given  values  of  speed  ratio,  the 
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following  relations  give  an  estimate  of  the  planar  surface 
of  solutions  that  appears  in  Figures  17  -  19. 

-  I'1  -  V  (610) 

where 

c  =  -89.40  degrees 
I  =  52.65  degrees 
a  =  -21.00  degrees 

and  values  of  a,  b,  and  c  are  polynomial  functions  of  u, 
given  by 

a  -  (106.701)  -  (271 . 537) u  +  (1209. 57)u2 

-  (1840. 30) u3  +  (895. 29) u4 

b  =  (5.657)  -  (90. 002) u  +  (401. 76) u2 

-  (61 1 . 66) u3  +  (297. 66) u4 

7o  =  (-26.104)  +  (81 . 784) u  -  (365.55)u2 
+  (556. 77) u3  -  (271. 00) u4 

Examples  of  planar  estimates  of  solutions  to  Eq  (6.5), 
generated  by  the  above  method,  are  presented  in  Figures  21  - 
23.  These  estimated  solutions  correspond  to  the  numerically 
generated  solutions  presented  in  Figures  17  -  19.  A  more 
exact  method  to  estimate  solutions  to  Eq  (6.5)  is  given  by  a 
three-dimensional  curvefit  algorithm  described  below. 


Figure  21.  Estimation  of  Trajectory  State  Solutions  of  Validity  for  the  Non 
Rotating  Earth,  dl/dh  Equation  of  Motion  (u  =  .860  ,  h  =  .0129) 
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Solutions . 


Using  polynomial  least  squares  curve  fitting 


routines,  a  curvefit  to  the  solutions  of  the  dl/dh  Rotate 
equation  was  developed  for  ranges  of  inclination  angle  and 
argument  of  latitude  at  epoch  of  0.5  I  <  75.0  degrees 
and  -45.0  a  <  0.0  degrees.  This  curvefit  is  three- 
dimensional  in  the  sense  that  it  predicts  solutions  for  7 
from  Eq  (6.5)  for  various  I,  a,  and  u.  Since  changes  in  h 
were  shown  to  produce  negligable  changes  in  7,  h  was  not 
included  in  the  fit.  The  resulting  algorithm  produces 
accurate  solution  estimates.  Figure  24  presents  numerically 
generated  values  of  7  versus  I,  and  Figure  25  presents 
corresponding  curvefit  generated  estimations  of  7  versus  I 
for  the  same  values  of  a,  u,  and  h. 
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Figure  25.  Estimation  of  Trajectory  State  Solutions  of  Validity  for  the  Non 
Rotating  Earth,  dl/dh  Equation  of  Motion  (u  =  .020  ,  h  =  .0052) 


VII.  Conclusions  and  Recommendations 


Conclusions 

This  effort  investigated  the  validity  of  the  non-rotating 
planet  assumption  for  three-dimensional  Earth  atmospheric 
entry.  This  study  was  limited  to  entry  at  orbital 
inclination  angles  between  0.5  and  75.0  degrees,  where  most 
Earth  atmospheric  entry  occurs,  and  vehicle  speeds  ranging 
from  circular  orbital  velocities  to  low  supersonic  speeds, 
where  terminal  maneuvers  (such  as  landing  approaches)  are 
usually  initiated.  Constant  lift-to-drag  ratio  and 
ballistic  coefficient  were  assumed  along  the  vehicle’s 
flight  path.  Validity  results  are  coordinate  dependent 
since  si ngu 1 ar i ti es  exist  in  the  equations  of  motion.  On 
the  basis  of  this  investigation,  the  following  conclusions 
are  made : 

1.  As  a  set,  the  five  non-rotating  planet  equations  of 

motion  (Section  IV)  are  invalid  for  Earth  atmospheric 
entry.  Hence,  the  non-rotating  Earth  assumption, 
common  in  analytical  entry  analyses,  produces 
incorrect  entry  trajectory  results.  Solutions  to  the 
non-rotating  Earth  equations  of  motion  were  derived 
from  expansions  of  the  rotating  Earth  equations  of 
motion  using  the  method  of  matched  asymptotic 
expansions.  These  solutions  are  given  by  Eqs  (5.86) 

-  (5.90) . 
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The  dq/dh  equation  of  motion  for  a  non-rotating  Earth 
is  never  valid  for  the  ranges  of  inclination  angle 
and  vehicle  speeds  mentioned  above.  Hence,  a 
solution,  Eq  (5.95),  was  developed  for  the  dq/dh 
equation  of  motion  for  three-dimensional  lifting 
atmospheric  entry  of  a  rotating  spherical  Earth. 

This  solution  is  valid  for  the  entire  non-dimensional 
altitude  domain  and  it  accounts  for  the  Coriolis 
acceleration  on  the  flight  vehicle.  The  method  of 
matched  asymptotic  expansions  was  employed  in  this 
derivation . 

A  variety  of  realistic  Earth  atmospheric  entry 
trajectory  states  exist  where  some  of  the  non¬ 
rotating  equations  of  motion  are  valid  for  a  rotating 
Earth.  The  dn/dh,  da/ dh ,  and  dl/dh  non-rotating 
equations  of  motion  are  valid  for  the  same  entry 
trajectory  states  for  a  rotating  Earth.  Other 
trajectory  states  exist  where  the  du/dh  non-rotating 
Earth  equation  of  motion  is  valid  for  a  rotating 
Earth.  These  two  sets  of  trajectory  states  do  not 
overlap  for  the  investigated  ranges  of  entry 
inclination  angle  and  velocity. 

A  number  of  trajectory  states  are  presented  where  the 
du/dh  non-rotating  equation  of  motion  and  the  dfi/dh, 
da/dh,  and  dl/dh  non-rotating  equations  of  motion  are 
valid.  Holding  speed  ratio  constant,  realistic 
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variations  of  non-dimensional  altitude  causes 
negligible  changes  in  these  states  of  validity.  The 
magnitudes  of  the  flight  path  angles  where  these  non- 
rotating  Earth  equations  of  motion  are  valid  increase 
as  speed  ratio  values  decrease.  Trajectory  states  of 
validity  are  symmetrical  about  various  values  of 
argument  of  latitude  at  epoch  for  the  non-rotating 
Earth  dl/dh  and  du/dh  equations  of  motion.  Two 
methods  to  estimate  flight  path  angle  values  where 
these  equations  of  motion  are  valid  are  also 
presented . 

Recommendations 

Based  on  the  observations  of  the  assumptions  and  of  the 
findings  of  this  investigation,  the  following 
recommendations  for  further  study  are  proposed: 

1 .  Investigation  of  the  validity  of  the  non-rotating 
planet  assumption  for  three-dimensional  Earth 
atmospheric  entry  was  limited  to  orbital  inclination 
angles  between  0.5  and  75.0  degrees.  Further 
investigation  should  be  undertaken  for  other  ranges 
of  inclination.  Near-polar  and  negative  inclination 
angles  would  be  of  particular  interest. 

2.  Further  study  of  the  trajectory  states  where  the 


du/dh  and  the  dfi/dh,  da/dh,  and  dl/dh  non-rotating 
Earth  equations  of  motion  are  valid  is  warranted.  An 


attempt  should  be  made  to  construct  entire 
trajectories  from  these  trajectory  states.  If 
successful,  a  reduction  in  the  complexity  of  the 
equations  of  motion  could  be  made  for  certain 
trajectories  and  possibly  even  specific  trajectory 
classes.  This  reduction  would  allow  the  non-rotating 
planet  equations  of  motion  to  be  used  in  place  of  the 
more  complex  rotating  planet  equations  of  motion, 
causing  past  analytical  studies  that  utilized  the 
non-rotating  planet  assumption  to  be  valid  for  these 
trajectories.  Current  trajectory  optimization 
computer  programs  consume  many  valuable  hours  of 
computer  run-time,  iterativly  integrating  the 
rotating  planet  equations  of  motion.  Optimization 
computer  programs  such  as  the  Air  Force  Flight 
Dynamics  Laboratory’s  ENTRAN  (ENtry  TRajectory 
ANalysis)  and  OTIS  (Optimal  Trajectories  via  Implicit 
Simulation)  codes,  could  greatly  benefit  from  a 
reduction  in  the  complexity  of  the  rotating  planet 
equations  of  motion. 

The  method  of  matched  asymptotic  expansions  has 
proven  to  be  a  powerful  tool  in  the  development  of 
solutions  to  boundary  layer  problems  such  as 
planetary  atmospheric  entry.  Further  application  of 
this  method  should  be  undertaken  to  obtain  a  full  set 
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of  solutions  to  the  equations  of  motion  for  three- 
dimensional  lifting  atmospheric  entry  for  a  rotating 
Earth.  This  activity  would  entail  finding  solutions 
to  the  complex  differential  equations  given  by  Eqs 
(5.20)  -  (5.24)  and  Eqs  (5.46)  -  (5.52),  and  then 
matching  and  forming  composite  solutions  from  these 


results . 


Appendix  A 


Derivation  of  Equations  Relating  6,  ,  j)  and  a,  0,  I 

Using  Spherical  Trigonometric  Relationships 


This  appendix  presents  the  basic  derivation  of  equations 
relating  6 ,  <p ,  ifi  and  a,  fi,  I  using  trigonometric 
relationships.  In  this  paper,  the  planetary  model  is 
spherical  and  hence  the  spherical  sine  and  cosine  laws  can 
be  successfully  applied. 

Figure  A1  depicts  a  reference  spherical  triangle.  For 
clarity,  the  angles  between  the  curved  line  segments  are 
called  interior  angles  whereas  the  angles  formed  between  the 


unit  vectors,  r.,  r .  ,  and  r,  ,  are  termed  exterior  angles. 

1  J  K 

Hence,  A,  B,  and  C  are  exterior  angles  and  a,  b,  and  c  are 
interior  angles  in  Figure  A1 . 

The  sine  law  for  spherical  triangles  states  that  the 
ratio  of  the  sine  of  an  interior  angle  to  the  sine  of  its 
opposing  exterior  angle  is  constant. 


sinCal  sin(b)  sin(c}  ,A  ,N 

sTn(A)  =  sin  (B)  =  iK(C)  (A1) 

The  cosine  law  for  a  spherical  triangle  states  that  the 
cosine  of  an  interior  angle  is  the  sum  of  two  products.  The 
first  product  is  formed  of  the  cosines  of  the  other  two 
interior  angles;  the  second  product  is  formed  of  the  sines 
of  the  other  two  interior  angles  and  the  cosine  of  the 


Figure  A1 .  Reference  Spherical  Triangle 
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opposing  exterior  angle.  The  cosine  law  is  given  by  the 
following  expression. 


cos  (a)  =  cos  (b)  cos  (c)  +  sin  (b)  sin  (c)  cos  (A) 


(A.  2) 


Eqs  (A. 1)  and  (A. 2)  are  used  to  relate  the  variables 
4>  >  0,  and  I.  From  comparison  of  Figures  11  and  12 

(Section  II)  to  Figure  A1 ,  the  exterior  and  interior  angl 
of  the  formed  spherical  triangle  are 


e, 


i 

i 

i 


es 


Exterior  Angles 

A  =  0  -  n 
B  =  4> 

C  =  a 


Interior  Angles 
a  =  k/2  -  f 
b  =  I 
c  =  k/2 
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where 


6 


longitude 

latitude 


<P  = 

ip  =  heading  angle 

a  =  argument  of  latitude  at  epoch 

0  =  longitude  of  the  ascending  node 

I  =  inclination  angle 

Application  of  the  sine  law  to  these  variables  gives  the 
following : 


i 


! 
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Si 
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s_in^  _  sing  _  sin(g  -  fl) 
sinl  sin(7r/2)  sin(7r/2  -  ip) 

Since 


cosip  =  sin(^  - 

*> 

and 

sin(|)  =  1 

sin(0  - 

sing  = - ^  — 

cos  ip 

HI 

(A.  3) 

simp  =  sinlsing 

(A.  4) 

Application  of  the  cosine  law  to  the  inclination  angle 
gives  the  following. 

cos  (I)  =  -cos  (^)  cos  -  r/>)  +  sin  (^)  sin  cos  (^) 

«■  «r 

Since  sin^  =  cos  (g  -  ip)  and  cos  (g)  =  0 

cos  (I)  =  cos(|t)cos(^)  (A.  5) 

Another  application  of  the  cosine  law  gives 

cos(^)  =  -cos(I)cos(^  -  ip)  +  sin(I)sin(^  -  ^)cos(g) 
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which  reduces  to 


cos(I)sin(^)  =  sin  (I)  cos  (y»)  cos  (a) 


or 

Yet 


cosa  = 

another  application  of 


tan^ 
tan  I 


the  cosine 


law  gives 


(A.  6) 


cos(^  -  i>)  =  -cos(|)  cos(I)  +  sin  (^)  sin  (I)  cos  (8  -  [}) 

sin(^)  =  sin(I)cos(0  -(])  (A. 7) 

Some  other  relationships  can  be  found  that  are  coupled 
to  Eqs  (A. 3)  -  (A. 7)  and  are  useful  for  some  occasions. 
Substituting  Eq  (A. 5)  into  Eq  (A. 6) 

cosa  =  cos^cos(0  -  n)  (A. 8) 

Substituting  Eq  (A. 5)  into  Eq  (A. 3)  gives 

sin<9  -  n>  =  IS!  <AB> 


Rewriting  Eq  (A. 6)  and  substituting  in  Eq  (A. 5)  gives  the 
expression 


sin  (ip) 


cosasinl 


cos  <p 

Substituting  Eq  (A. 4) 


into  this  result 


The  relations  given  by  Eqs  (A. 3)  -  (A. 10)  are 
Section  III  to  transform  the  equations  of  motion  f 
of  latitude,  longitude,  and  heading  angle  to  terms 
containing  orbital  inclination  angle,  longitude  of 
ascending  node,  and  argument  of  latitude  at  epoch. 


(A. 10) 

used  in 
rom  terms 

the 
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sin(eb)  - 


Therefore 


3,3  5,5  7,7 

e  b  ,  e  b  e  b 

3!  5!  ~  7!  + 


sin(eb)  =  eb  +  0(e  ; 

As  an  example,  the  following  expansion  of  the  cosine 

2 

function  is  derived  to  Q(e  ). 

cos  a  +  eb  +  e2cj  =  cos  (a)  cos  |eb+e2  cj  -  sin  (a)  sin  Jeb  +  e 

=  cos (a) ^cos (eb) cos (e2 c)  -  sin (eb) sin  (e2  c)  J 

2  2  1 

-  sin(a)  sin(eb)cos(e  c)  +  cos(eb)sin(e  c) 

o 

=  cos  (a)  -  eb-sin(a)  +  0(e  ) 

The  following  relations  were  obtained  through 
application  of  the  Taylor  series  expansion: 

cos  j^a  +  eb  +  0(e2)]  =  cos(a)  -  eb’sin(a)  +  0(e2) 

sin  a  +  eb  +  0(e2)j  =  sin(a)  +  eb’cos(a)  +  0(e2) 

tan  a  +  eb  +  0(e2)l  =  tan(a)  +  — — -  +  0(e2) 

J  cos2  (a) 


tan [a  4  eb  +  0(e2)]]  =  ^ 


tan  (a)  .  2  ,  >. 

'  sin  (a) 


4 -  +  0(£S) 


[ein[*  -  eb  *  0(e!)]]  =  sdrno  -  £b’— 1 Hfr  4  0(€l) 

'  sm  (a) 


*.•  ^  "  - 
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Appendix  C 


Derivation  of  Solutions  to  Selected  Outer  Expansion 


Differential  Equations 


Outer  Expansion  £  Terms 


The  e  terms  of  the  outer  expansions  found  in  Section  V 


are  as  follows: 


dh  "  (1+h) 


(C.l) 


%  % 

(1+h) L  u 


(C.2) 


(C.3) 


(C.4) 


dh  (l+h)tan7 


(C.5) 


Solutions  to  this  set  of  differential  equations  are 


derived  below. 


du/dh  Equation  Solution.  Rearranging  Eq  (C.l)  gives 


o  -dh 


(1+h) 


or  ln(u  )  =  -ln(l+h)  +  K 


Solution : 


u  (1+h)  =  C 
0  v  '  1 


(C.6) 


_%  V**  "»  "■  _%  .*1  — *»  %  **»  ~w  j,”**  — ^  fc  >  ■  I,  *  a  •>  | 
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dq/dh  Equation  Solution.  From  Eq  (C.2) 


dqo  %  \  %  ,1 

dh  _  (1+h)  u  ~ 


Application  of  Bernoulli’s  equation  to  this 
differential  equation  makes  the  solution  derivation 
straightforward.  Bernoulli’s  equation  is  given  by 

^  +  P(x)y  =  Q(x)yn 


This  form  gives  the  solution  (Beyer,  1984:315) 

r 

i/exp  [  (l-n)  Jp(x)dxj  =  (1-n)  Q(x)  [  (1-n)  Jp(x)dxj  dx  +  C 


where  v  =  y  n  and  n  ^  1 


lienee  ,  +  P(h)%  =  Q  C^)  % 


where 


n  =  3,  v  =  q^ 


P(h)  =  iTThy  ’  and 


ech)  >  h-  =  nW 

i  o  v  y 


Solution:  JyT+hy  =  ^-n(^+^) 

r 

q^  2  exp  -21n(l+h;j  =  ^  exp  j^-21n  (1+h)  j  dh  +  K 

1  - 


(l+h)l 


f 

•  (1+h) _2  dh  + 


Integrating  gives 


% 


V  =  §  (1+h)  -  c2(l+h)2 


(C. 


or 


%  =  [§  (1+h)  -  C2  (1+h)2] 


-1  /  2 


(C. 


dl 

_ c 

dh 


dnc 

dh 


dl/dh 

Equation 

Sol ution . 

From  Eq 

(C.3) 

=  0  ; 

I 

=  C 

0 

5 

dn/dh 

Equation 

Solution . 

From  Eq 

(C.4) 

=  0  ; 

n 

=  C 

0 

A 

da/dh 

Equation 

Solution . 

From  Eq 

(C.5) 

(C. 


(C. 


da. 


dh  (l+h)tan7 


By  definition,  q  =  cos' y  .  The  expansions  of  q  and  C0S7 
given  by  the  following  expressions: 

qQ+  eq]  +  0 (e2 )  =  cos  7q  +  £7^  0 (e2 ) ] 


=  cos70  -  e7jsin70+  °(e  ) 


By  noting  the  order  of  e  in  this  equation,  it  can  be  seen 
that 


7) 


8) 


9) 


10) 


are 


%  =  cos70 


(C.ll) 


From  Eqs  (C.ll)  and  (C.12),  an  expression  is  found  for  the 
tangent  of  7  . 


tan  7  =  — — - 

'0  COS7 


(1  -  cos!7o)’/2  (1-q,2)1/2 


tanTo  =  [“T  ”  *] 


(C . 13) 


lere  was  previously  found  to  be 


=  J-(l+h)  -  C2  (l4h)2] 


-1  /  2 


(C . 14) 


Therefore 


k  -  c*(i+h)2  -  i]"'/2-(ttet 


Making  x  =  (1+h)  causes  dx  =  dh  and  allows  this  integral 
to  be  rewritten  as 


%  =  j[§.  'x  -  V2  - 


-1  /  2 


The  integral  is  now  in  the  form  of 


-  dx 

°°  =  Jx(X),/2 


(C • 15) 


where  X  =  a  +  bx  +  cx2  ,  f  -  4ac  -  b2  .  k  =  4c/f . 


■vv 


i  • 


f  =  2^C\  C2  ~  1)'  and  k  =  “°1  C2/(C1  °2  "  1} 

c 

1 


The  solution  to  Eq  (C.15)  is 


ao  = 


dh 


x  (X) 


1  /  2 


=  -  (-a)'l/2cos  !  [  (hx  +  2a)/{lxl  (-q)1 /2>]  +  < 


Hence 


cos 


-l 


2x/C  -  2 


21x1  (C2-l/C  2)1/2 


+  C. 


This  solution  can  be  expressed  in  either  of  the  following 
two  f orms : 


a  =  -  cos 


-l 


1  -  C  /(1+h) 


(1  -  c  2c2)l/2 


+  c„ 


(C.ll 


or 


u  =  1  +  cos 


(a  -  C  )  •  f  1  -  C  2  C  1 
v  0  3  J  l  1  2  J 


1  /  2 


(c.r 


Outer  Expansion  e  Terms 


The  e1  terms  of  the  outer  expansions  found  in  Section 
IV  are  as  follows: 


du  -u 


l  1 


dh  ”  (1+h) 


]  1  /2 

3u_ (1+h) I  • cosl 


0 


0 


(C  •  1 


dl 


dh 


=  2 


3 (1+h) 


u 


1  /  2 


cosa  sinl 


~  ~ - -(cosa  tan7  -  sina  )  (C.l 

tan7  v  0  '0  o ' 
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dir  +  TTTKy  =  -4cos(C5).(3C 

This  differential  equation  is  in  the  form  of 

* 

f 

^  4  P  (x)  y  =  Q  (x) 

This  form  gives  the  solution 
y '  exP  [/p(x)dx]  =  Jq(x)  [/p(x)dx]dx  4  C 
For  the  du^ /dh  differential  equation 
y  =  u  >  X  =  h, 

P(h)  =  ~(T+h)  »  and  Q  =  -  4cos(Cb  )  *  (3C]  ) 

exp[/(f;hy]  =  exp[ln(1+h)]  =  (1+h) 

(l4h)U]  =  qJ  (l4h)dh  =  |(l+h)2  4  C6 

Sol ution : 

C  / 

u,  =  O^h)  ~  2(3C)1/2cos(Cb)(l4h) 


dq/dh  Equation  Solution. 


where  uQ  ,  u^  ,  and  1^  are  given  by  the  following  equations: 


C  /(1+h) 


(C.6) 


-  2(3C)1/2cos(C5)(l+h) 


(C . 23) 


I  =  C 
o  b 


(C.9) 


Substitution  of  Eqs  (C.6),  (C.9),  and  (C.23),  into  Eq  (C.22) 


gives 


dq,  V 

dh  ~  Z  C 


%  qi 


2  (1+h) 


-  2q()  (l+h)cosC&  [|-] 


+  (3C  ),/2  (l+h)2cosC 

C  1  1  b 


(C . 27) 


iere  qQ  =  §(l+h)  -  C?  (l+-h)2] 


-1  /  2 


(C.8) 


Eq  (C.8)  can  be  substituted  into  Eq  (C.27)  to  produce  a 
complicated  integral  with  h.  However,  for  typical  lifting 
earth  atmospheric  entry,  altitudes  of  interest  range  from  a 
minimum  of  0  to  a  maximum  of  about  100  nautical  miles.  This 
is  equivalent  to  a  range  of  non-dimensional  altitude  of 
about  0 . 0  h  <  0.029  .  The  maximum  h  encountered  for  a 

typical  lifting  entry  trajectory  for  Earth  is  thus  always 
less  than  the  small  parameter,  e,  allowing  h  to  also  be 
utilized  as  a  small  parameter  in  the  outer  expansion 
differential  equations.  The  use  of  h  as  a  small  parameter 


KWV 


&E5S5 


v.  rjsrTJsiKv  ’wv»  rj>  tj<  r>  'ji  -*j<  ■•,'• ', 


considerably  simplifies  the  derivation  of  solutions  to  some 
of  the  first  order  e  differential  equations. 

Expanding  the  terms  in  Eq  (C.27)  for  small  h  and 

neglecting  the  resultant  terms  of  0(h2)  greatly  simplifies 

the  analytical  integration  of  this  equation.  One  term  that 

appears  often  in  Eq  (C.27)  is  qQ .  Rewriting  Eq  (C . 8) 


a  =  [p-  +  ^  -  C  -  h2C  -  2C  hi 
o  Lt  t  2  2  2  J 


-1  /  2 


1  1 


A  condition  allowing  the  expansion  of  this  equation  for 
small  h  is  given  by 


§“(l+h)  >  C2  (1+h)2 


This  condition  can  be  easily  seen  to  always  hold  true  for 
any  values  of  C^  ,  ,  and  h. 

Expanding  Eq  (C.8)  for  small  h  using  a  Taylors  series 
approximation  gives 

%  =  [e;  -  CJ  ’/!-  h[e7  -  c2] ' [§7  -  202]  +  0(h2)  <c-28) 

Some  other  terms  appearing  in  Eq  (C.27)  are  likewise 
expanded : 


su(1+h)  =  1  + 


o  _  n  n 


one 


r 

+  Oi 


(C .  29) 


V 


=  [§--c2J  -  3h[§-  -  c2]/  IS-  -  C2J  +0(h2)  (c-31) 


Defining  two  new  constants,  both  dependent  upon  and  C^  , 


simplifies  Eqs  (C.28)  -  (C.31). 


let  C,  J  =  [§ - C2  ]  and  C,  ,  =  [t-  -  C,  ] 


With  substitution,  Eqs  (C.28)  -  (C.31)  become 


%  =  (C,  ,) 


-1  /  2 


C  h 

1  2  +  0  (h2  ) 


(C  ) 

v  1  l  J 


l  /2 


(l+h)  =  (C  ,)'l/2  + 


l  1 


c,  <c, , > 


3  /2 


+  0(h2) 


,  1  2C1!h  2 
V  =  ^ - -  0(h  ) 


CM  (C,,)’ 


V  -  <c,,)'3/2-  7?TT72  +  0(h2) 


(C  ) 
v  1 1  ' 


Two  other  terms  in  Eq  (C.27)  can  also  be  expanded 


V 


(1+h)2  =  (Cn)  3/2  + 


(C  ) 
v  1  1  J 


3  /2 


3C 


2  - 


1  2 


1  1 


[TT  =  -q.  +  q.h  +  0  (h2  ) 


(C . 32) 


(C . 33) 


(C .  34) 


(C . 35) 


(C . 36) 


+  0(h2)  (C .  37) 


(C .  38) 


Substitution  of  Eqs  (C.33)  -  (C.38)  into  Eq  (C.27)  gives  a 


solvable  expression  for  dq  /dh . 


<*q 


L  =  q,  C,  a  +  qCr  h  +  C  h  +  C 


dh  ‘‘l  16  15”  14”  13 


(C .  39) 


where 


16  C  (0  )'/2  C,Cu 

1  v  1 1 J 


-  1 


(C .  40) 


2C 


ci6  “  1  C  C 


1  2 


1  1  1 


1  1 

Cn  +  (C  1)l/2 


(C . 41) 


3C  C 
6  1  2 


14  (C,)2(Cn)5/2 


+  2cos (C  ) 

D 


1-3 


1  2 


1  1 


c, 3  c, , 3 


1  /  2 


(C .  42) 


C,  3  =  2cos(C5) 


1 

3 

c  c  1 

C  C 

1  1 1  J 

l  l  l  J 

T  1  /2 


3C 


C  2C  3/2 

l  i  l 


(C . 43) 


The  differential  equation  given  by  Eq  (C.39)  is  in  the  form 

+  P(x)y  =  Q(x) 

This  form  gives  the  following  solution 
y  ’  exp  |/p (x) dxj  =  Jq(x) [/p(x) dxj dx  +  C 

For  the  dq^ /dh  differential  equation  y  =  q  ,  x  =  h, 


With  substitution 


exp[Jp(h)dh]  =  exp[-Cigh  -  §(c15)h2]  =  exp[_C16h  4  °(h2)] 


Therefore 


r  r 

qiexp[~Ci6h]  =  c!3  exp[-C16h]dh  +  C  4  h[exp(-C  6h)]dh  +  C 


and  the  solution  to  the  dq  /dh  equation  is  given  by 


qi  "  C 


1  3 


1  4 


-  h 


1  6 


(C  )• 
V  1  6  ' 


1  4 


LC1  6  J 


-  C  exp[-C  hi 

20  KL  ) 6  J 


(C. 


The  constants  C  ,  C  ,  and  C  are  given  by  Eqs  (C . 43) , 
1314  16 


(C.42),  and  (C.40),  respectively.  Because  li  terms  are 


neglected,  C  is  dropped  from  the  solution.  Hence,  a 


l  b 


solution  for  q  has  been  found  as  a  function  of  two 


constants  of  integration,  Cg  and  C^Q  ,  and  the  constants  C 


C  ,  and  C  ,  which  are  dependent  upon  C  and  C  from  the 
14  3  6  12 

relations  given  in  Eq  (C.32). 
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Appendix  D 

Derivation  of  Solutions  to  Selected  Inner  Expansion 


Differential  Equations 


Inner  Expansion  eu  Terms 


The  e°  terms  of  the  inner  expansions  found  in  Section  V 
are  as  follows: 


o  d  1  Lcosa 

0  sin70  cdcos70 


(D.l) 


^ 


Be  *  cosff 


(D.2) 


sinacosa 


cos7osin7o 


(D.3) 


sinasina 

_ o 

cos7q  sin7o  sinIQ 


(D.4) 


smasma 


tanl  sin^  cosl 

0  '0  0 


(D .  5) 


Solutions  to  this  set  of  differential  equations  are  derived 


below . 


dq/d£  Equation  Solution.  From  Eq  (D.2) 


L  n  -f 

-  ■  Be  'cosa 

CD 


This  is  a  very  simple  integration  by  separation  of 
variables.  The  constant  of  integration  is  K  . 


(E 


a  =  rf— -  Be  Vosa  +  K 
^0  cD  1 


du/dh  Equation  Solution.  From  Eq  (D.l) 


du 

37^-  =  -  2u  Be  ^ 
d£  0 


1  L . COSCT 


sin70  cd  cos7 


0  J 


Previously  it  was  seen  that  q^  =  cos7o 


T.  ,  %  d1o  .  d1o  CL„  -f 

Therefore,  ^7—  =  -sin).—  and  3 =  7,— Be  s 


0  d£ 


dS  "  CD 


The  du/dh  equation  can  be  rewritten 


^“o  - 


-2u  Be 


0  sm) 


0 


-2u  Be  f 


0 


7)  % 


dJ±.d£_ 

d?  d')0 

de 


fi.df, 
d«  d% 


ln(u0)  = 


1  +1 


-2u  Be 

0  Sin7 


i  -1 


L.Be-£  cosa 
CD  sm7 


i  +1 


-2u  Be“« 

0  CD  %  J 


CL  „  -5 

—  Be  '•cosa 
CD 


0  J 

-1 


d7r 


lqo 


After  reducing,  this  equation  takes  the  form 


ln(u0)  = 


f-CP2dt  + 

Cj^COSCT 


2  d 
% 


a  =  <y2l°  +  21n(7  )  + 

T)  C^COSCT  'O 


Solution : 


u0  =  K3%  exp 


■CD27o 


Cj^cosct 


dl/d£  and  da/d£  Equations .  F rom  Eq  (D . 3) 


COSCT 

sin70 


In  (K3  > 

(D 


a.L,  oT  .  sinacosa 

—2.  =  _k.Be_£ - °- 

di  CD  cos70sin70 


From  Eq  (D.5) 


CT  *  sinasina 

-  p— 'Be  £  •  —  — - — — j 

G~.  tanl  sinnr  cos! 

D  0  '0  0 


By  the  chain  rule 


dl  dl  d£ 

da,  df  da, 

0  s  0 


tanl 
_ 0 

tana. 


With  separation  of  variables,  this  integration  simply 
produces 


sina  sinl  =  sinK 

0  0  4 


s  i  nK  "* 

I  =  sin"1  — - -I 

0  sino. 


da/d£  and  dn/d£  Equations.  From  Eq  (D.5) 


^  sinasma0 

tanl  sinnr  cosl 

0  '0  0 


From  Eq  (D.4) 


dn  CT  f.  sinasina, 

o  =  _L.R  -C - - - -°_ 


dC  ‘  % 


cos7q  sin7o  sinIQ 


Chain  rule 


daQ  daQ  d£ 

dn  d£  dn. 


=  -cosl 


Substitution  gives 


da,  ,  rsinK  ,1  rsinK  ,2' 

jn  =  -cos  sin  [ — : - 1  =  -  1  -  [ — t - 1 

dn,  lsina„J  lsinaAJ 

0  .  0  .  0 


Rewriting  this  equation 


sina  da 

0  0 


an  = 


[sin2a0-  8in2Kj 


l  /  2 


To  solve  this  integral,  let  x  =  cosa  and 


a  =  1  -  sin2K  =  cos2K  , 
A  A 


0  9 

so  dx  =  -  sina  ’da  and  1  -  x  =  sin  a 

0  0  0 


Therefore,  the  integral  can  be  written  in  the  form 


dx 


J  [a2  -  x2]1/2 


=  -  cos 


■i  r_x_i 

llal 


Solution:  [}Q  =  -  cos  -  K, 


or 


cosa  =  cosK  cos (K  -  0  ) 

U  i  o  0 


dl/d^  and  d^y/d^  Equations.  From  Eq  (D.3) 


dl  CT  ..  cosa„ 

_ 0  q  c  •  0 

— —  =  *=— •  Be  ssina- 


d£ 


"D 


cos' v  sin^ 


o 


From  above 


d7n 


d£ 


CT  f 

L.Be-C  cosa_ 


Bin7, 


o 


By  the  chain  rule 


dIo  dIo  d*  .  cosao 
d7n  dC  d7n  cos7n 


Since 


.  -i  rsinK<i 

a,  =  sin  v--T~ 


then  by  substitution 


tana 


cos7„ 


rsinK  , ' 

cos  sin  1  — : — ?4-  dl 

Lsinl  J  0 

J  L  0  J  J 


This  can  be  rewritten  as 


tana 


COS7 


sin2  K.  1  1/2 


•  2  T 
sm  X. 


sinl 

_ 0 

.  [sin2  -1  +  s 


hsin  I 


2  T  1 1  /2  0 


To  solve  this  integral,  let  x  =  cosl  and 


2  i  •  2  t/  2  r/ 

a  =  1  -  sin  K  =  cos  K 
4  4 


so  dx  =  -  sinl  dl  and 
0  0 


1  -  x2  =  sin2I„ 


Thus,  the  right  side  of  the  integral  can  be  written  in  the 


and  the  left  side  of  the  integral  is 

*  j 

tan<T  =  Analog  [tan  [|  +  J-]  ]  +  Kg 

0 


Solution : 


cosIQ  =  cosK^ cos [tana  *  log [tan + 
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Inner  Expansion  £  lerms 

The  e1  terms  of  the  outer  expansions  found  in  Section  V 
are  as  follows: 


CT  -  f-  a ,  sina, 

p— -Be  ^sina  - - - - 

CD  [cos70sin7o 


tana  dl 
o  .  i 


'1  2 
COS  7 
»n 


rl  •  2 
sm 


tanl  d£ 


dOj  s 

ar  =  crBe  tsl™  % — ? 

*  U  cos 


(D . 13) 


(D . 14) 


sma„  sma 

0  0 

- 2 -  -  7j  ~  2 - 

cos  7  tanl  sin  7  tanl 

*00  '00 


sma„  cosa„ 

-  _ 0 _  +  a  _ 0 _ 

"1  .  2  T  costl  sin7  tanl 

cos7o  sin7o  sin  IQ  '0  '0  0 


o  „  -£  cos70  cosff 

=  2u»Be  ■»,— — +  eraser 

sm  D  '0 

•n 


!t  Ct  7  cosa 

L  cosa  L  '1 


Cr>  •  2 

D  sin7  cos  7 


2u^  Be 


1  Lcosq 

sin70  +  Cdcos70 


(D.15) 


dn  tana  dl 

_ L  _  _ 0  .  1 

d£  sinIQ  d£ 


dni  CL  „  -f  .  [  sinao 

ar  =  cr  Bs  sin<T  t — 5 — rr 

s  D  cos  7^  sinl 

•  O  1 


'  1  •  2 

sinl  sin  7  sinl 

0  0  '0  0 


smasina,  cosa, 

_  1  _ 0 _  +  7, _ 2 _ 

1  .  .  2-r  COS7,  sin7  sinl 

cos70  sin7o  sm  1^  '0  '0  0 


(D .  16) 


(D. 17) 


i|'  |4» 
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P 

e 

s 


'2?: 


To  solve  the  dq  /dh  equation  is  a  trivial  exercise; 


solutions  to  Eqs  (D.13)  -  (D.16)  are  not  required. 

dq/d^  Equation  Solution.  From  Eq  (D.17),  q^  is  simply 


q,  =  K? 


(D.18) 


This  appendix  presents  three  example  Earth  atmospheric 
entry  trajectories  to  give  an  indication  of  what  the 
dimensionless  variables  u  and  h  are  in  relation  to  more 
conventional  variables  such  as  Mach,  altitude,  flight  path 
angle,  and  velocity.  These  trajectories  were  generated  by 
AFWAL/FIMG,  Air  Force  Flight  Dynamics  Laboratory  personnel 
in  the  course  of  flight  performance  analysis  support  for 
current  hypersonic  entry  vehicle  studies.  Modern  trajectory 
analysis  computer  programs  were  utilized  to  produce  this 
data . 

Example  Trajectory  1  is  the  gliding  entry  trajectory 
for  a  relatively  high  lift-to-drag  ratio  vehicle  (about 
3.0).  Example  Trajectory  2  is  the  gliding  entry  trajectory 
for  a  low  lift-to-drag  ratio  vehicle  (about  0.5).  Example 
Trajectory  3  is  a  very  steep  gliding  entry  trajectory  for 
the  same  low  lift-to-drag  ratio  vehicle.  This  last 
trajectory  is  somewhat  non-realistic  due  to  the  extremely 
high  dynamic  pressures  (and  hence  severe  heating) 
encountered.  However,  it  is  included  to  help  give  an 
indication  of  the  extreme  in  u  and  h  combinations.  Many 
current  and  planned  lifting  entry  vehicles  have  trajectories 
that  fall  within  the  ranges  of  u  and  h  shown  in  these 
examples . 


For  the  sample  trajectories  presented  in  Tables  I,  II, 
and  III,  the  following  variables  are  defined. 

y  =  Altitude  (kft) 

7  =  Flight  path  angle  (deg) 

V  =  Velocity  (kft/s) 
h  =  Non-dimensional  altitude  =  r/r 


u  =  Speed  ratio 
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262.3 

227.9 
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138.9 
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107.0 


Table  II  Example  Traiector 
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Table  III  Example  Trajec 
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Figure  F6.  States  of  Validity  for  the  Non-Rotating  Barth 
du/dh  Equation  of  Motion  (u  =  .210  ,  h  =  .00' 
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Figure  F9 .  States  of  Validity  for  the  Non-Rotating  Earth 
dl/dh  Equation  of  Motion  (u  =  .210  ,  h  =  .00 


Figure  F10.  States  of  Validity  for  the 
dl/dh  Equation  of  Motion 


States  of  Validity  for  the  Non-Rotating  Earth 
dl/dh  Equation  of  Motion  (u  =  .010  ,  h  =  .00 


Figure  F12.  States  of  Validity  for  the  Non-Rotating  Earth, 

dl/dh  Equation  of  Motion  (u  =  .210  ,  h  =  .0077) 


Figure  F13.  States  of  Validity  for  the  Non-Rotating  Earth, 
dl/dh  Equation  of  Motion  (u  =  .210  ,  h  =  .007 
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The  assumption  of  a  non-rotating  planet,  common  in  most  analytical  entry 
trajectory  analyses,  has  been  shown  to  produce  significant  errors  in  some  solutions 
for  the  lifting  atmospheric  entry  of^Earth.  This  thesis  presents  an  investigation 
70^' the  validity  of  the  non-rotating  planet  assumption  for  general  three-dimensional 
Earth  atmospheric  entry. 

In  this  effort,  the  three-dimensional  equations  of  motion  for  lifting  atmospheric 
are  expanded  to  include  a  rotating  planet  model.  A  strictly  exponential  atmosphere, 
rotating  at  the  same  rate  as  the  planet,  is  assumed  with  density  as  a  function  of 
radial  distance  from  the  planet's  surface.  Solutions  are  developed  for  the  non¬ 
rotating  Earth  equations  of  motion  and  for  one  of  the  rotating  Earth  equations  of 
motion  using  the  method  of  matched  asymptotic  expansions. 

It  is  shown  that  the  non-rotating  Earth  assumption  produces  incorrect  entry 
trajectory  results  for  entry  orbital  inclination  angles  between  0.5  and  75.0^  > 
degrees  and  vehicle  speeds  ranging  from  circular  orbital  velocities  to  low 
supersonic  speeds.  However,  a  variety  of  realistic  trajectory  states  exist  where 
some  of  the  non-rotating  Earth  equations  of  motion  are  valid.  Three  of  the  non¬ 
rotating  equations  of  motion  are  found  to  be  valid  for  the  same  entry  trajectory 
states.  Other,  independent  trajectory  states  exist  where  a  fourth  non-rotating 
Earth  equation  of  motion  is  valid.  A  fifth  equation  of  motion  is  never  valid  for 
the  ranges  of  orbital  inclination  angle  and  speeds  investigated Trends  in  the 
results  of  the  trajectory  states  of  validity  are  discussed  and  niAthods  to  estimate 
some  of  these  states  are  presented.  / 
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